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THEORETICAL STUDY OF CONDITIONS LIMITING V/STOL TESTING 
IN WIND TUNNELS WITH SOLID FLOOR 

By Harry H. Heyson 
Langley Research Center 

SUMMARY 


Theoretical study indicates that, under sufficiently severe conditions, the forward 
portion of the wake will flow forward along the floor. The location of the consequent 
rolled-up vorticity results in the testing limits found previously by Rae (Journal of 
Aircraft, May- June 1967). The predominant factor determining the initial flow reversal 
of the wake is wake deflection. Separation of the flow on the ceiling of a closed tunnel is 
also indicated, but at substantially more severe conditions. While it is doubtful that 
wind-tunnel data corresponding to free -air conditions can be corrected successfully 
beyond Rae's limits, it may be possible to obtain wind-tunnel data for ground-effect 
conditions well beyond these limits. 


INTRODU CTION 

Many of the theoretical assumptions which result in comparatively simple correc- 
tions for wind-tunnel tests of conventional aircraft are inadequate for tests of V/STOL 
models. The principal assumptions of classical theory that are violated in V/STOL 
testing are that the wake is undeflected from the free-stream direction and that the wall- 
induced velocities are only small perturbations on the free-stream velocity. 

To a certain extent, the aforementioned assumptions have been removed from more 
recent theoretical treatments (refs. 1 to 4) intended specifically for use with V/STOL 
models. On the other hand, recent experimental evidence (particularly ref. 5) indicates 
that there are definite limits beyond which the theory of references 1 to 4 cannot be 
applied. 

In an effort to explore the reason for these limits, the present investigation exam- 
ines the entire flow in the wind tunnel rather than only that portion induced by the tunnel 
boundaries. The analysis assumes a simple lifting rotor as representative of a very 
high-lift, large-downwash lifting system. For computational simplicity, the rotor wake 
is presumed to pass downward and rearward in a straight line until it meets the floor, 
after which the wake is assumed to flow off along the floor. Even though such an ideal- 
ized wake could not exist in its own induced field, it is shown that the assumed wake does 
indicate, at least qualitatively, the correct experimentally known wake deformations. 


Thus, examination of the predicted deformations in the tunnel yields qualitative informa- 
tion regarding the effects of the walls on the general flow within the tunnel boundaries. 

It will be shown that the limits proposed by reference 5 appear to be associated 
with a region of reversed flow on the floor of the tunnel. The effects of a number of 
parameters (wake skew angle, wind-tunnel width-height ratio, height of the model above 
the floor, ratio of model size to tunnel size, and angle of attack (representing changes in 
drag-lift ratio)) on the region of reversed flow are studied in some detail. 

Certain preliminary results of this investigation, for one rotor and one wind-tunnel 
configuration, have been published previously as reference 6. 

SYMBOLS 


The symbols used in this paper largely conform to those of references 1 to 4 with 
only minor exceptions. Extensive use is made of direct output from digital computers 
and digitally controlled plotters. The available symbols and lack of lower case letters 
in the output of these devices necessitate deviations from the symbols used in the text. 
The alternate versions, indicated in parentheses, show the forms used in machine- 
prepared tables and figures. The appendixes of this paper use symbols which may differ 
from those of the main text. Separate symbol lists are provided, where required, in the 
appendixes. 


At 

B 

b 

C D (CD) 
C L (CL) 

C Lhb 

C-O-B 

D 


momentum area of lifting system 
cross-sectional area of test section 
semiwidth of test section 
complete span of model 
generally, induced drag coefficient, D/qS 
lift coefficient, L/qS 

lift coefficient based on wind-tunnel area below model and within span, 
L/qhb 

abbreviation for closed-on-bottom-only 
generally, induced drag 
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H 


semiheight of wind tunnel 


h height of model above floor 

L lift 

1 2 

q (Q) dynamic pressure, — PV 

Cl 

q (QC) dynamic pressure corrected for boundary interference 
R rotor radius 

S area of lifting system 

u,v,w induced velocity components, parallel to X-, Y-, and Z-axes 

V forward velocity or uncorrected wind-tunnel velocity 

V c wind-tunnel velocity corrected for boundary interference 

Wq mean, or momentum theory, value of induced velocity parallel to Z-axis, 

positive when directed upward 

X,Y,Z in wind tunnel : Cartesian coordinate system centered in wind tunnel at 

model; X positive rearward; Y positive to right; Z positive upward 
for rotor : Cartesian coordinate system centered in, and oriented with 
respect to, rotor tip-path plane; directions as above 

x,y,z (X,Y,Z) distances from origin along X-, Y-, and Z-axes 

x^ distance behind model of the intersection of the wake and the floor 

a (ALPHA) rotor tip -path-plane angle of attack, positive when upward in front 

a c corrected angle of attack 

y (GAMMA) width-height ratio of wind tunnel, B/H 

Au boundary-induced interference velocity parallel to X-axis, positive rearward 
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Aw boundary-induced interference velocity parallel to Z-axis, positive upward 

Aa (D ALPHA) boundary-induced additive correction to angle of attack 

£ (ZETA) ratio of tunnel semiheight to the height of the model above floor, 

H/h 

7 ? (ETA) in references 2 to 4, distance from model to right-hand tunnel 

boundary divided by tunnel semiwidth (Throughout the present 
paper, ri = 1.) 

p mass density of test medium 

a (SIGMA) ratio of rotor radius to tunnel semiwidth, R/B 

X (CHI) wake skew angle, angle measured from negative Z-axis to center 

of wake, positive rearward 

RESULTS AND DISCUSSION 

DISTRIBUTION OF WALL-INDUCED INTERFERENCE 

In general, wind-tunnel-interference calculations omit the direct effect of the model 
and are concerned only with the additional effects created by the near proximity of the 
wind-tunnel boundaries. At the outset it will be instructive to examine the manner in 
which the boundary-induced interference velocities vary throughout the wind tunnel. 

In the present study, which assumes that the model is a simple lifting rotor, the 
interference velocities are easily obtained using the computer program given as appen- 
dix K of reference 4. This program provides the interference at a tail behind a rotor; 
simple modifications (to set the tail span equal to zero and to enclose the tail length and 
height within DO-loops scanning the X-Z plane of the tunnel) produce the results shown 
in figure 1 for a series of decreasing wake skew angles. In figure 1, the rotor has uni- 
form load distribution, is mounted 30 percent of a semiheight above the center of a 
square, closed test section, with a = 0.6 and a = 0°. For each wake angle, contours 
of equal values of the, ratios Aw/wq and Au/wq have been obtained by noting that, in 
references 1 to 3, A m /Ax = ^ c fiy. 
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Interference at x = 90° 

Figure 1(a) shows the interferences at x = 90°, that is, when the wake passes 
directly rearward as in classical wind-tunnel-interference theory. Observe that wq is 
always negative for positive lift. Thus, except for a few regions far forward of the rotor 
(where the rotor produces upwash in free air) the vertical interference is always in the 
nature of an upwash which increases with distance downstream. Since the rotor is 
mounted above the centerline, the interference is greater above the rotor than below the 
rotor. (A centered rotor would produce a field symmetric about the X-Y plane.) The 
longitudinal interference is significant only near the model and results from the effective 
bound vortex represented by the rotor. An increase in velocity above the rotor and a 
decrease below the rotor are evident, with the values above the rotor being somewhat 
greater than those below the rotor. (This pattern would be antisymmetric about the X-Y 
plane if the rotor were centrally located.) 

Now the effect of these interference velocities is to alter the effective angle of 
attack and forward velocity of the rotor; thus, from reference 2 


where 


and 


so that 


Oq = a + Aal 


A a = 



1 


+ 


Au 

V J 



( 1 ) 


( 2 ) 


Examination of equations (1) and (2) indicates that it is the ratios with respect to V 
rather than the ratios with respect to Wq which are of primary concern. The ratios 
with respect to V are obtained by observing that 


Aw _ /Aw \ / w o\ 
v >o/W 



(3) 
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where, for the conditions of figure 1 (a = 0°), equation (A8) of appendix A shows that 

w n 

~ = -cot x (4) 

Strict application of equations (1) to (4) indicates that the effect of the interference 
at x = 9°° is identically zero. This result is expected since, with Wq = 0 (eq. (4)), 
there is no lift on the rotor. 

The result at x = 90° is not, however, trivial. Provided that Wq is small and 
that the wake does not depart by more than a few degrees from the horizontal, the actual 
value of w 0 /v, rather than the theoretical value (eq. (4)), must be used in equations (1) 
to (3) to obtain usable results. It will be observed that this treatment of wall effects for 
the limiting case of x = 90° is the procedure used in most wind-tunnel-interference 
theory. 


Interference at x = 70° 

Similar interference-velocity contours for a more heavily loaded operating condi- 
tion (x = 70°) are shown in figure 1(b). It may be noted that the wall -induced interference 
grows rapidly as the distance along the wake increases. At the intersection of the wake 
and the floor, it is obvious that the vertical interference must be exactly equal, but oppo- 
site to, the vertical velocity produced by the rotor in free air, otherwise there would be 
flow through the floor. The longitudinal interference at the rotor remains small, but this 
interference also increases rapidly with distance along the wake. At the floor, the longi- 
tudinal interference opposes the tunnel velocity ahead of the wake and augments it behind 
the wake. Following equations (3) and (4), the magnitude of the longitudinal interference 
is found to be about 11 percent of the tunnel velocity for the largest contour (Au j wq = ±0.3^ 
shown in figure 1(b). 


Interference at x = 50° 

As the rotor lift is increased to produce a wake skew angle of 50° (fig. 1(c)) the 
primary effect on Aw/wq is merely an increased gradient along the wake. This is the 
expected result since the effect of the floor is still to just cancel the free-air field, but 
the length of the wake within the tunnel is reduced. A similar increased gradient is 
noticed for the longitudinal interference; however, in this case the levels of interference 
at the floor also increase to approximately twice those obtained at p 70°. Together 
with the increase in w qJv, the net result is longitudinal interferences at the floor which 
are about 50 percent as great as the tunnel velocity. 
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Interference at x = 30° and 10° 

Figures 1(d) and 1(e) show continuations of the same trends as the wake steepens. 
The maximum indicated longitudinal interference velocities on the floor become about 1.4 
and 5.0 times the tunnel velocity for x = 30° and 10°, respectively. 

In the face of such extreme wall -induced interferences, it is obvious that the local 
flow (and the wake, which follows the local flow) near the floor must undergo drastic 
deformations. The magnitude of the anticipated deformations is so great that quantita- 
tively accurate theoretical results can not be obtained by the present methods; however, 
qualitatively, it appears that the theory does indicate the correct trends to be observed 
in the flow patterns. The remainder of this paper will be concerned with an examination 
of the nature and effects of these deformations. 

INDUCED FIELD OF A LIFTING ROTOR 

The complete flow within the wind tunnel consists not only of the boundary-induced 
interference velocities, but also the general tunnel velocity and the velocities induced 
directly by the model. The tunnel velocity presents no difficulty since, in the theory, it 
merely requires the addition of a constant term to the calculated flow. The induced field 
of the rotor is more involved, and its calculation will be examined in this section of the 
paper. 


Method of Calculation 

For many purposes, it is adequate to represent a rotor wake as a skewed ellipti- 
cally cylindrical surface of uniformly distributed vorticity. This wake model has been 
used extensively in the past (refs. 7 to 10) and has been shown experimentally (refs. 9, 

11, and 12) to yield usable results for the time-averaged induced field of the rotor. Even 
in ground effect, where the real wake may be extremely distorted when compared to the 
simple cylinder, usable qualitative results are obtained (ref. 13). Since the present study 
is qualitative in nature, the simple cylindrical wake is used herein, thus avoiding the 
computational complexity inherent in a more detailed representation. 

The component of induced velocity normal to the rotor disk (w) is available in the 
form of nondimensional charts in several papers (for example, refs. 8, 9, and 14). No 
similar charts are available for the u and v components of the field. In any event, 
for the present purpose, where the field at large distances from the rotor is required, 
it is necessary to obtain the values by the direct numerical integration of equations (9), 
(11), and (13) of reference 10 (with i(\p) in that paper set equal to 1.0 so as to obtain 
the field for a uniformly loaded disk). 
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Typical Calculated Field 


A sample of the field calculated in this manner is presented for the X-Z plane of a 
rotor (x = 50°) in figure 2. Along the X-axis, the normal component (w) of induced 
velocity changes from an upwash ahead of the rotor to an increasingly large downwash at 
the trailing edge of the rotor, after which the downwash decreases asymptotically to zero 
at an infinite distance downstream. Passing downward along the wake, it is observed that 
the downwash increases and essentially achieves its ultimate uniform value of w/wq = 2 
within a diameter of the rotor. The longitudinal component (u) is of the nature of a sink 
in the region above the rotor and draws the flow into the wake. Within the wake, the flow 
is directed rearward, essentially achieving its ultimate value of u/wq = 2 tan(x/2) 

(refs. 6 and 7) within a diameter of the rotor; however, outside the wake and below the 
rotor, the flow is generally retarded. In particular, in the lower left quadrant of figure 2, 
the longitudinal component of rotor-induced velocity is such as to augment the previously 
discussed tendency of the wall-induced velocities to oppose the tunnel velocity. 

Effect of Disk- Load Distribution 

In general, the actual rotor disk-load distribution is far from uniform (ref. 9) and 
the nonuniformity of loading has strong local effects on the distribution of induced veloc- 
ities in the rotor field. This effect will be examined briefly in a later section of this 
paper, but the computed results presented herein are largely for uniform loading. Justi- 
fication for this usage is obtained from reference 11 which indicates that the field of a 
uniformly loaded rotor is a close approximation to the experimentally determined aver- 
age field over a span equal to that of the rotor. 

COMPLETE FLOW IN WIND TUNNEL 

Calculated Flow Based on Vortex-Cylinder Wakes 

Calculation procedures .- The complete flow in the wind tunnel is the sum of the 
wall-induced interference, the rotor-induced field, and the free-stream velocity. The 
accuracy of the calculation, at least for the assumed wake model, is dependent upon the 
consistency with which the image wakes used to represent the boundaries represent the 
wake within the tunnel. It is desirable to use the same representation for all of the 
wakes, both interior and exterior to the test section. 

Appendix B presents the derivation of an extended version of reference 1, which 
maintains the desired consistency of representation but which, in its original form, suffers 
from several defects. In common with reference 1, the theory presented in appendix B 
is restricted to zero angle of attack since this restriction is inherent to the methods of 
superposition used in the derivation. The longitudinal and lateral interference velocities 
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are incorporated in the revised theory by an application of the equations of reference 10. 
Appendix B also revises the image system of reference 1 so that the rotor may assume 
any arbitrary vertical position in the tunnel. Appendix C presents the FORTRAN pro- 
gram used to obtain the results presented herein. 

Calcul ated results .- The theoretical results obtained herein are presented in the 
form of vectorial flow charts prepared on a digitally controlled plotter (CALCOMP). 
Appendix D presents the most general of the FORTRAN programs used to prepare con- 
trol tapes for the automatic plotter. This program is illustrative of the general 
technique. 

Vectorial flow charts of the complete flow are presented in figures 3 to 30 for the 
same conditions as in figure 1. In addition, similar charts are presented in figures 3 
to 30 for the equivalent conditions in free air, in ground effect, and in a wind tunnel 
closed only on the bottom. Included in figures 3 to 30 are the flows in the X-Z plane, in 
the plane of the floor, and in a transverse plane through the center of the intersection of 
the wake and the floor. Isometric projections of the flow are also provided. Table I 
summarizes the conditions considered and provides, in addition, the calculated correc- 
tions, both to free air and to ground effect, as obtained from references 2 to 4. 

Before examining the computed flows it should be noted that reference 15 (see also 
refs. 6 and 9) has shown that the effect of wake rollup is to effectively modify the wake 
skew angle to some greater mean value. This feature has not been included in the pres- 
ent analysis since it leads to certain inconsistencies in depicting the entire flow. The 
present study is intended to be qualitative only and the effect of wake rollup will be mini- 
mized herein by comparing flows corresponding generally to the same wake angles. 

The general growth of the rotor free-air flow field with decreasing skew angle 
(increasing Cl) is evident in figures 3 to 30. Except at the lowest wake angle (10°) 
there is an upwash at the leading edge of the rotor disk. This implies that the initial 
path of the front part of the wake will be upward, returning back downward through the 
disk shortly thereafter. The flow at the trailing edge of the disk is sharply downward, 
and for high lift coefficients (y < 40°) is actually directed forward to some extent as well. 
In the longitudinal (X-Z) plane the net result is contraction of the actual wake. These 
features are actually observed in flow visualization studies. (See, for example, ref. 16 
and fig. 19 of ref. 2.) 

The effect of imposing additional restraints on the flow depends, of course, on the 
nature and location of the restraint. Thus, the rotor -induced vertical velocities are can- 
celled at the ceiling of the closed tunnel and augmented at the upper boundary of the 
closed-on-bottom-only tunnel. At an equivalent height above the rotor, ground effect is 
hardly noticeable. It is noteworthy that the flow at the ceiling can actually be reversed 
for sufficiently severe conditions. (See fig. 9(c) in which y = 10°. ) In practice, such a 
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flow reversal would lead to large-scale separation of the flow from the ceiling, with con- 
sequent extreme effects on the entire flow (ref. 17). The large downwashes for similar 
conditions at the upper boundary of the closed-on-bottom-only tunnel should lead to sim- 
ilar concern about the actual behavior of the main wind-tunnel stream (parts (d) of 
figs. 14 to 16). 

The floor is solid for all three constraint configurations considered herein, and the 
effects are strikingly similar as a result. Within the primary rotor wake the flow gen- 
erally turns sharply rearward at the floor; however, as the wake steepens below angles 
of 30° an increasing portion of the forward wake turns forward. This behavior is 
required in order to approach the symmetrical conditions which must exist in hovering 
(y = 0°). Large lateral velocities are observed at the sides of the wake near the floor 
(figs. 10 to 30) indicating that the wake will spread laterally at the floor. 

Behind the wake near the floor, there are a general acceleration of the flow and an 
upwash that are in sharp contrast to the general deceleration and downwash that are pres- 
ent in free air. The calculated results indicate that the actual rear portion of the wake 
probably will, in reality, turn rearward well above the floor and be carried off in such a 
manner as to approach the floor only asymptotically, if at all. 

The flow near the floor in front of the wake is the most complex and the most 
interesting. At very mild conditions, such as x - 70° (figs* 3, 17, and 24), the flow in 
this region, although slightly slower than in free air, is still predominantly rearward. 
Thus, the forward portion of the wake will turn rearward at the floor, and the entire wake 
will travel off above the floor. Under such conditions, the wake will, in all probability, 
never actually reach the floor. Instead, it will sink until the self-induced vertical veloc- 
ities of the wake are in equilibrium with the boundary-induced upwash and thereafter will 
continue in distorted form but with no further downward progress. This behavior is veri- 
fied by experiment (ref. 16). 

When the wake skew angle is depressed to 60° (figs. 4, 18, and 25), the calculated 
flow ahead of the wake at the floor is observed to be essentially stagnant. The forward 
edge of the wake, which is carried off with the average of the velocities just inside and 
just outside the vortex sheet, will still travel rearward; however, the turn rearward would 
be expected to be sharper and closer to the floor. With more sharply inclined wakes 
(figs. 6 to 9, 19 to 23, and 26 to 30), the calculated flow in this location is actually directed 
forward, counter to the prevailing tunnel velocity. The strength and extent of the reversed 
flow region grows as the wake angle steepens. In addition, the rearward flow just within 
the wake edge weakens. Finally, at some wake angle in the vicinity of y = 40° (figs. 6, 
20, and 27) the average longitudinal velocity across the forward edge of the wake passes 
through zero. For sharper wake angles, below 30° (figs. 7 to 9, 21 to 23, and 28 to 30), 
the presence of a large standing vortex on the floor is indicated. 


10 



On the floor, except near the side boundaries, the flow is essentially identical 
(figs. 17 to 23) for all constraint configurations considered herein. This result occurs 
because, at such close proximity, the image directly below the floor, which represents 
ground effect, is responsible for such large interference velocities that the sum of the 
effects of all the other images (provided to represent the top and sides of the tunnels) is 
insignificant by comparison. Near the side constraints, however, there are significant 
differences (figs. 10 to 30). Compared to ground effect, the flow in the closed tunnel is 
squeezed inward, and the flow in the closed-on-bottom-only tunnel is drawn outward. 

Effe ct of wake distortions .- The foregoing analysis, which is predicated upon the 
assumption of a rigid wake, has indicated the presence of a large degree of wake distor- 
tion. These distortions must now be examined to determine whether or not the changes 
in the field caused by the distortions augment or diminish the trends indicated by the sim- 
plified theory. 

Figure 31 presents a sketch of the deformed and undeformed wakes near the floor. 
The direction of vorticity in the real wakes and their images in the floor is indicated. In 
the rear (and also near the sides of the wake for which a similar sketch could be drawn) 
it is evident that the velocities will be generally increased at the floor and that the 
increased velocities will cover a larger region of the floor than for the undisturbed case. 

In the forward region the flow reversal is augmented. The extent to which the dis- 
tortion adds to the computed velocities as the wake changes from a rearward deformation 
at low lift to a forward deformation (as indicated in fig. 31) should be particularly large. 
As a result, the actual reversal in direction of the wake as Cl increases should occur 
with a "snap" action. Considerable hysteresis might be expected; that is, if the reversal 
resulted from a small increase in Cl> a substantially larger decrease in Cl might be 
required to remove the reversal. 

In the presence of a thick boundary layer, considerable local unsteadiness might be 
expected at the change in direction of the wake. This unsteadiness may be only local and 
not be of any consequence at the model. Both before and after the wake reverses there 
will be a steady flow. Thus, there may be no unsteadiness noted at the model at all, and 
the presence or absence of unsteadiness is no criterion as to the presence or absence of 
the phenomenon. 

Most of the features of the flow that have just been discussed have been observed 
experimentally. (See, for example, refs. 18 to 20.) 

Fin al form of wake .- The forward-going portion of the wake is overrun by the free- 
stream velocity and is immediately below a region of upwash. Consequently, a large 
standing vortex will form on the floor in front of the wake. Even the present rigid-wake 
calculations indicate this behavior (figs. 7 to 9 and 28 to 30), and the effect of wake dis- 
tortion is to magnify the result. 
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The general nature of the flow on the floor (figs. 17 to 30) is such as to transport 
this rolled-up vorticity sidewards and rearward along a more-or-less parabolic path. 
This path will be somewhat narrower in the closed tunnel, and somewhat wider in the 
closed-on-bottom-only tunnel, than in ground effect. In ground effect, the path of the 
vorticity will continue smoothly beyond the region for which results are presented in fig- 
ures 3 to 30. In the tunnel with open sides (parts (d) of figs. 3 to 30), the consequent path 
of the vorticity after reaching the sides of the tunnel is somewhat problematical and this 
path can probably be determined only by experiment. In the closed tunnel the vorticity, 
upon nearing the walls, will be predominantly influenced by its own image in the walls 
(not present in the theoretical treatment) which will drive the rolled-up vorticity up the 
walls as it passes rearward. A sketch of the resulting wake pattern is presented in fig- 
ure 32. 

Many experimental studies indicate that the indicated wake shape is correct. The 
flow studies of reference 16 clearly show the rolled-up vorticity flowing upward along the 
sides of the tunnel; unfortunately however, the windows of the tunnel used for the studies 
reported in reference 16 were too small to show the behavior on the floor. Reference 19 
presents sketches of the observed flow, which, in the fully developed state, are essentially 
identical to figure 32. (See, in particular, fig. 10 of ref. 19.) 

Reference 21 presents tuft-grid flow studies of a helicopter rotor in a 30 by 60 foot 
(9 by 18 meter) closed-on-bottom-only wind tunnel. Figure 33 presents a comparison of 
the calculated flow and the observed flow. The general nature of the field is represented 
reasonably well by the theory, with the experimental rollup, as expected, being somewhat 
more severe than indicated by the theory. At the location of the tuft-grid, about 2 feet 
(0.6 meter) to the side of the center of the tunnel, the flow should be essentially similar 
to that in either ground effect or in a closed tunnel. 

Although the treatment herein is specifically for rotors, it is noted that the same 
general results should be obtained for any high-lift system. This observation is verified 
by wind-tunnel studies of fan-in-wing models (ref. 19), rotors and lifting propellers 
(refs. 5 and 20), and jet-flap models (refs. 15 and 18), as well as ground-effect studies of 
tilt-wing (ref. 22) and deflected- slipstream (ref. 23) aircraft. In addition, similar results 
have been observed in unpublished studies by Eldridge (see ref. 24) utilizing jet-lift 
models. 

Effect on wi nd-tunnel data.- A wind-tunnel flow pattern, similar to that of figure 32, 
is obviously greatly different than the flow in free air. Consequently, the forces mea- 
sured on the model will not be representative of those which would be obtained in free 
air, and it is not possible to directly relate the measured data in any simple manner to 
an equivalent free-air flight condition. 
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Since data taken under conditions of this nature are neither valid nor truly useful, 
it is obvious that there is some limit beyond which tests of a particular model in a par- 
ticular tunnel cannot be obtained. This limit does not necessarily correspond to the 
initial behavior of the leading edge of the wake; for example, the large lateral velocities 
at the edge of the wake near the floor could, for narrow tunnels and wide models, con- 
ceivably start the rolled-up vortices of figure 32 even before the leading edge of the wake 
snaps forward. 

The effect at the model of the rolled-up vortices will depend upon their strength 
(determined by some combination of model size, model forces, and tunnel velocity) and 
their location with respect to the model (determined by some combination of model size, 
tunnel size, tunnel velocity, and tunnel proportions). One parameter which indirectly 
incorporates many of these test conditions has been shown in reference 24 to correlate 
the "limit” data of reference 5. This is the ratio of the impingement distance of the wake 
on the floor behind the model to the model span, or 

(5) 

The actual correlation is shown in figure 34 where it may be seen that x-fjb, for the 
data of reference 5, falls into two bands (dependent upon the presence or absence of large 
fillets in the tunnel corners) which are a function of the wind-tunnel proportions only. 

Two additional points indicated on the figure show that essentially the same limiting val- 
ues are obtained for a jet-flap model (ref. 15) and a tilt-wing model (ref. 25). 

The data of reference 5 were all obtained near zero drag-lift ratios (a = 3°, so that 
D/L ~ 0.05), and little experimental evidence is available at other drag-lift ratios. The 
point shown for the jet-flap model was obtained at D/L ~ 0.3 and provides an indication 
that, at least, in this range D/L has little effect other than its influence on y. (See 
eq. (A5) of appendix A.) A certain amount of caution should be used in extrapolating to 
extremely different conditions in using equation (5). 

Equation (5) indicates that testing may be extended to higher lifts or lower speeds 
if the height of the model above the floor is increased. Experimental confirmation of this 
result has been obtained recently (ref. 20). 

Corrections for e xtreme conditions .- Conceptually at least, it is possible to include 
the actual self-induced wake distortions in the wall-interference theory and to obtain 
the wall-induced interference as the difference between the free-air and wind-tunnel 
fields. Initial efforts in this direction for moderately heavily loaded wings have already 
been published (refs. 26 and 27). Under conditions as severe as those considered herein, 
the analysis would be considerably more involved than in recent so-called free-wake 
analyses (such as ref. 28). Even presuming that convergent numerical results were 
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obtained, there is considerable doubt that the large field distortions obtained could be 
converted into satisfactory data corrections (refs. 5 and 24). In any event, such an effort 
would appear uneconomic in relation to the storage and speed available in current digital 
computers. 


Calculations Based on Doublet-String Wakes 

At this point it would be desirable to survey a number of conditions other than the 
one treated so far. Unfortunately, such a survey is not economically feasible with the 
program of appendix C which, on the CDC 6600 computer, requires approximately 
0.3 minute per point (or 5 hours for fig. 24). A faster, although also somewhat limited 
method, is developed and used in the succeeding portions of this paper. 

Calculation procedures.- Fundamentally, reference 2 develops wall effects by con- 
sidering the limiting case of a vanishingly small model as represented by a wake that is 
an inclined line density distribution of doublets. Finite-size models can be treated in this 
method by superposing the computer interference fields of a suitable assemblage of such 
doublet-string wakes (refs. 2 to 4). 

Figure 35 shows the pattern of doublet-string wakes used in references 3 and 4 to 
represent a rotor of finite size and compares this assemblage to the uniform skewed vor- 
tex cylinder employed in appendix B. In the limit, if an infinite number of doublet strings 
distributed uniformly through the cylinder were employed, the field would be identically 
that of the cylinder. Indeed, a more elaborate version of this technique has been used to 
calculate rotor blade loads (ref. 29). However, when only 20 strings are employed, as in 
references 3 and 4, the field within and very near the wake will be somewhat "lumpy” 
compared to that of the vortex cylinder. Normally, however, distances this close to the 
wake are not involved in wall-effect calculations. 

The treatment of references 3 and 4 is not limited to zero angle of attack as is the 
theory presented in appendix B, Although identical superposition techniques are used 
(see ref. 2), the individual semi -infinite doublet strings, each representing only a small 
part of the wake, are independently truncated at the floor. Nonuniformity of load distri- 
bution requires no significant increase in computer time since the effect of a nonuniform 
load distribution is represented by a simple rearrangement of the assemblage of doublet 
strings (ref. 3). 

Because of the lumpiness of the field of the doublet-string model, the velocities 
produced by the rotor itself must be obtained from the vortex- cylinder model. For angles 
of attack other than zero, the wake skew angle, distances, and velocities, all of which are 
determined from an axis system embedded in the tip-path plane, must be properly 
adjusted to correspond with the wind-tunnel axes. 
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Since the representation of the real rotor wake is different from that of its 
images, erroneous results will be obtained at the immediate intersection of the real wake 
and the image wake below the floor. Since none of the doublet strings in the arrangement 
used in references 3 and 4 lie on the longitudinal axis of the wake, the computed velocities 
on this axis are relatively lower than they should be. Thus the combination of the two 
representations on the floor should indicate a flow through the floor. This effect will 
cover the greatest region when the wake and its image graze along the floor at high wake 
skew angles and, the affected region will diminish with decreasing skew angles. In the 
flow charts presented herein the most obviously incorrect of the flow vectors affected by 
this mismatch have been deleted, although vestiges of the aforementioned behavior 
remain. 

Because of the deleted vectors near the floor there is a substantial loss of precision 
in the most interesting region of the flow. Certain conclusions can be drawn, however, 
from the nearby flow. 

The theoretical treatment of references 2 to 4 does not include the lateral compo- 
nent of wall-induced interference. While this component could be obtained by differenti- 
ating equations (1) and (7) of reference 2 with respect to y and then proceeding as in 
appendix B, the effort does not appear to be productive because of the aforementioned 
imprecision in the most interesting region of the flow. Instead, the present results will 
be confined to the X-Z plane of the tunnel where symmetry requires that the lateral 
interference velocities are zero. 

The required modifications to the program of appendix K of reference 4 are 
straightforward with the real wake being incorporated from reference 10 by means of 
subroutines RTRL0AD and ROTVEL which are given as part of appendix C in this paper. 
Consequently, a complete listing of the program is not included herein. The complete 
program, running on the CDC 6600 computer requires about 1/2 second per point (com- 
pared to 0.3 minute for the program of appendix B). Results for a triangular disk-load 
distribution require about 10 percent more time, solely spent in evaluating the field of 
the real wake. 

Figures 36 to 125 present flow charts, based on a doublet-string representation of 
the wake, for a wide variety of operating conditions and rotor and wind-tunnel configura- 
tions. Table EL, which serves as an index to these figures, also presents the calculated 
average boundary-induced interference as obtained from appendix H of reference 4. 

Comparison of results from t wo methods .- Figures 37, 40, 43, and 45, prepared 
using the doublet-string wake, are otherwise identical to figures 3, 5, 7, and 9, respec- 
tively, which were prepared using vortex-cylinder wakes. Comparison of the corre- 
sponding figures of these sets indicates that the results are practically identical. 
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Effect of constraint configuration. - Although there are significant differences 
throughout the remainder of the flow, the flow near the floor is essentially identical for 
the three constraint configurations. (Compare parts (b), (c), and (d) of each figure.) The 
reason for this behavior is the consistently closed lower boundary in all three 
configurations. 

Effect of wake skew angle.- Comparisons at equal angles of attack, between fig- 
ures 36 to 125 indicate that the dominant parameter affecting the flow reversal on the 
floor ahead of the wake is the wake skew angle. Irrespective of the values used for any 
of the other parameters in the calculation, this flow is never reversed for a skew angle 
of 70° and is always reversed for a skew angle of 30°. With some allowance for the dif- 
ferent location of the constant grid of vector locations with respect to the wake under 
varying conditions, it would appear that the flow is also always reversed at a wake skew 
angle of 50°. 

The actual rollup of the rotor wake is ignored in the calculated flow. Since the 
actual wake follows a rearward curving path as it rolls up and since this path should, for 
an otherwise constant model, be a function of only the model dimensions, it would be 
anticipated that the actual local wake skew angle near the floor will be larger for small 
models than for large models. Further, since the predominant influence on the flow is 
the image in the floor, and even more specifically, those portions of this image nearest 
the floor, a certain interrelationship between size and wake angle should be observed in 
the real wind tunnel. Although the wake angle will still be the predominant influence, the 
required wake skew angle, as computed at the rotor, may be somewhat smaller for 
smaller models. 

Effect o f angle of attack (D/L).- The effects of large excursions in angle of attack 
(corresponding to large changes in D/L) are shown in figures 36 to 125 to be relatively 
small for a constant wake angle. This observation is particularly true when the large 
requisite changes in force coefficients are considered (table n). At constant wake angle, 
the strength and extent of the region of reversed velocity ahead of the wake are slightly 
decreased for negative angle of attack and slightly increased for positive angle of attack. 

On the other hand, the rearward velocities within the wake are increased substan- 
tially at negative angles and correspondingly decreased at positive angles. Since the point 
at which the wake flips forward rather than rearward depends upon the average longitudi- 
nal velocity on each side of the vortex sheet, a somewhat steeper wake angle will be 
required for wake reversal at large negative ou Consequently, at wake reversal, the 
reversed region ahead of the wake may, by virtue of the steepened wake angle, be larger 
and stronger than the corresponding flow at a = 0°. 

The behavior to be expected of the real curved wake is unclear. The vastly 
increased force coefficients at negative angle of attack will result in a more rigid wake 
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with a larger angle of inclination with respect to the floor. This effect would tend to 
counter the effect of rollup. The actual behavior probably can only be determined from 
well-planned experiments. 

Effect of model vertical location .- Figures 36 to 65 illustrate the flow for three 
different vertical locations of the model above the wind-tunnel floor. There is no effect 
evident on the initial appearance of reversed flow on the floor. The vertical extent of the 
reversed flow region diminishes softiewhat with decreasing model height; however, the 
decrease in model height is more than adequate to maintain the rotor within the region of 
maximum disturbance at the lower wake angles. 

Considerations of rollup, however, lead to the conclusion that the local wake angle 
at the floor will be smaller at the higher model locations. Thus, when real wakes are 
considered, some delay in flow reversal may be expected for high model mountings. This 
result is indicated from the experiments of reference 20; however, the results presented 
therein probably are influenced to a greater degree by the subsequent path of the rolled- 
up vorticity than they are by the initiation of wake reversal. 

The flow at the ceiling of the closed tunnel is affected greatly by model height. The 
tendency of the flow to reverse, and consequently to separate, from the ceiling is 
increased substantially by raising the model. Thus, there are limits to which the effects 
of flow reversal can be modified by changes in model height. Too high a model height 
may produce more severe effects in the tunnel than those which the experimenter is 
attempting to avoid. 

For the range of model heights considered herein, ceiling separation always 
requires more severe testing conditions than those required to initiate flow reversals. 
Thus, within these limits, tests conducted with larger values of xf/b than those pre- 
sented in figure 34 will most likely be safe from the occurrence of ceiling separation. 

Effe ct o f rot or load distri bution .- The flows shown in figures 66 to 75 for a trian- 
gular load distribution may be compared directly with figures 46 to 55 to determine the 
effect of load distribution. Although major differences due to loading are evident within 
the wake, any differences in the region of the floor in front of the wake are insignificant. 

Effe ct of relati ve r otor size .- Figures 76 to 85, 46 to 55, and 86 to 95 present 
sequences of flow charts which differ only by uniform increments in the ratio of rotor 
diameter to tunnel width. The overall interference level in the tunnel increases almost, 
but not quite (see ref. 3), as rapidly as the ratio of rotor area to tunnel area (ita^y/A). 

The reversed flow region at the floor is noticeably stronger and more extensive for the 
larger rotors; however, there is no significant effect on the wake angle at which the cal- 
culated flow initially reverses. As noted earlier, there will be some interaction between 
rotor size and the local wake angle at the floor (due to rollup) which will result in a 
somewhat earlier appearance of reversal for larger rotors. 
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Effect of w ind-tunnel prop ortions .- Figures 96 to 115 present flow charts for a 
width-height ratio (y) of 2.0 with a equal to 0.6 and 0.3. Figures 116 to 125 present 
similar charts for y = 0.5 and <7=0.6. Comparisons between these charts and fig- 
ures 36 to 95 indicate that the net effect of width-height ratio on the flow at the floor is 
negligible. The effects shown in figures 96 to 126 are almost purely those which would 
have been expected on the basis of the interrelated effects of y and a in altering 
Ajh/At (=7rcr2y/4) and R/h (=oy£). As noted earlier, changes in tunnel proportions will 
alter the eventual path of the rolled-up vorticity and, in consequence, change the effect on 
the data. 


COMPARISON WITH REFERENCE 18 


Results of Reference 18 

Using tufts and smoke streamers, reference 18 has studied the conditions under 
which ’’flow separation” initially occurs on the floor of a closed wind tunnel for jet-flap 
models. The correlation of such conditions is presented in terms of a lift coefficient 
Cuhb b ase d on the tunnel cross-sectional area below and within the wing span. The 
results are a function of the model drag-lift ratio, where, being the total balance force, 
the drag is the total drag. (Note that there will be some difference expected between this 
drag-lift ratio and the drag-lift ratio, defined in terms of induced drag, used in the 
remainder of the paper. This difference is ignored herein, since the profile portion of 
the experimental drag, including separation, cannot be estimated adequately.) The 
resulting critical values of Cl^ are given in reference 18 as 


c Lhb = 3 *° 



(6a) 




(6b) 


Techniques of Reference 18 

Since reference 18 is concerned solely with the appearance of anomalous effects on 
the tunnel floor, it would appear at first glance to be ideal for correlation with the pres- 
ent theoretical study. Unfortunately, several features of the techniques used in refer- 
ence 18 tend to cloud the comparison. 

First, the results of reference 18 are based on tufts on the floor and smoke 
streamers introduced along the floor. Thus, at best, the point of "initial separation" 
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most likely coincides with the initial appearance of a reversed flow region ahead of the 
wake rather than with the actual reversal of the wake itself. Information as to the 
behavior of the wake is absent, although it could have been obtained by introducing smoke 
into the wake at the model. 

Secondly, there was a boundary layer of unknown thickness and condition present on 
the floor during the tests. As shown in reference 30, it is possible for the disturbance at 
the floor to propagate forward through the low-energy flow of the boundary layer with 
marked effects on the resultant flow. One solution to this problem is to eliminate the 
boundary layer on the flow by a technique such as that of references 31 and 32; however, 
such an approach would have been clearly beyond the scope of reference 18. 

Finally, the models, which were all jet-flapped wings were somewhat limited in 
performance. Large negative drag-lift ratios can only be obtained at very low tunnel 
velocities and with large negative angles of attack. The low speeds present considerable 
difficulty in measuring the actual tunnel velocity, particularly in the presence of tunnel 
flows which may be reversed locally; the large negative angles of attack imply leading 
edge separation with a consequent "bluff-body 1 ’ wake and large blockage. 

Nevertheless, no alternate study is available which addresses itself so closely to 
the effects treated herein. The foregoing effects should be kept in mind throughout the 
subsequent discussion. 


Comparison at D/L = 0 

For zero drag, equation (A21) of appendix A shows that may be expressed 

as 


^ _ 7r(b/h) 

c Lhb — L - 


tan x sin x 


(7) 


Equation (7) may be solved for b/h with Clj^ = 3.0 to obtain the limiting span- 
height ratio at zero drag as 


U Q 

U = i ^ x sin x 


( 8 ) 


Equation (5), used earlier to correlate the data given by Rae in reference 5, simi- 
larly may be solved for the maximum allowable span-height ratio as 


b _ tan x 
h ( x f/ b ) 


(9) 
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If it is noted that Xf/b ranges between about 1.25 and 2.0 (fig. 34), and that the 
wake angles corresponding to Rae's limits (ref. 5) are on the order of 30° to 60°, equa- 
tions (8) and (9) appear to lead to the same order of maximum allowable span. The com- 
parison, however, is probably not legitimate since South in reference 18 was concerned 
solely with an initial disturbance at the floor, whereas Rae in reference 5 was concerned 
with determining the point at which wind-tunnel disturbances had grown so large that his 
force data was affected to a noticeable degree. The disparity of definition is obvious in 
equations (8) and (9). In agreement with the present theoretical study, equation (8) shows 
that there is no dependence whatever on the wind-tunnel proportions. As discussed pre- 
viously, the effect on the data will be a function of the path of the entire rolled-up vortex 
system within the tunnel and this path in turn will be a function of the wind-tunnel propor- 
tions. This dependence is inherent in equation (9) since Xf/b is a function of the tunnel 
width -height ratio (fig. 34). 

South's experiment (ref. 18) is more nearly comparable with the present study. 
Observe that for zero drag, equations (7) and (8) indicate that South's limit is primarily 
dependent on wake skew angle. A dependence on the span is also noted. As discussed 
previously (under the heading "Effect of Wake Skew Angle"), a similar dependence would 
be expected from the theory if wake rollup were considered. 

Comparison at Positive D/L 

Figure 126 compares South's limit (eq. (6b)) with the calculated values of C^b 
(eq. (A21)) for a series of constant wake skew angles and three values of b/h. Irre- 
spective of b/h, it may be seen that, for positive D/L, South's limit corresponds to a 
constant wake skew angle within the accuracy of the data of reference 18. (The data 
given in that paper scatter about the correlating curve by about ±0.5 in terms of C]_,hb') 
There is, of course, a relationship between the wake angle and the span-height ratio in 
figure 126. At least in direction, this relationship is that which was anticipated earlier 
by considering the effect of rollup on the local wake angle near the floor. 

Comparison at Negative D/L 

The constant limiting value of C^^ = 3.0 for negative D/L is in glaring con- 
trast to the behavior of the limiting values obtained under positive D/L conditions. 

There is no dependence of C^b on the wa k e skew angle, whereas, even considering 
wake rollup, the theoretical study presented herein indicates that the flow reversal should 
still occur at a nearly constant wake angle. 

The following analysis will indicate that the limits proposed in reference 18 can not 
be applied to conditions well beyond those examined in that paper; that is, to conditions 
representing, for example, a propeller in nearly axial flight. 
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Figure 127 shows the theoretical behavior of (for b/h = 1.0) at negative 

drag-lift ratios as great in magnitude as 100. Observe that for any wake angle other 
than 90°, C L hb always becomes infinite at some value of D/L. Equation (A21) of 
appendix A yields this result since C^b will approach infinity as the term tan x + D/L 

approaches zero. Since tan x = - — - — (eq. (A5) of appendix A), it is obvious that the 

w 0 L 

infinite C^b occurs when 


0 = tan x + — 
L 


_V_ _ D + D _ __V_ 
Wq L L Wq 


or 


V = 0 


when 


tan x = 


D 

L 


( 10 ) 


This result would be anticipated since Cj^^, being based on dynamic pressure, must 
become infinite under static thrust conditions. Equation (10) merely states the static 
thrust conditions in terms of x and D/L. Since C^t) approaches infinity for 
D/L = -tan x, it must certainly pass through any finite value first. The required value 
of D/L at constant x and C^b is obtained by solving equation (A21) for D/L to 
yield 


D/L = -tan x ± \j(V C Lhb)( b /h)sec x (11) 

where the positive sign corresponds to forward flight and the negative sign corresponds 
to rearward flight which is not under consideration herein. 

Figure 128 compares the values of D/L for zero velocity and for C^b = 3.0. 

As anticipated, the constant value of CLhb (obtained for less negative D/L in ref. 18) 
would imply a flow reversal or "separation" on the floor of the tunnel for all static thrust 
conditions and for some speed range above V = 0, irrespective of model size, even if the 
model (a propeller, for example) were only slightly misalined with respect to the flow. 

It should be observed that reference 18 notes that 'The initial stages of separation 
did not necessarily produce a significant flow up the tunnel side walls in the vicinity of the 
model" (as in fig. 32). Subsequent discussion between the present author and the author 
of reference 18 has brought out that the quoted behavior was generally the case for neg- 
ative D/L, but that, in general, abrupt and violent flow up the walls (as in fig. 32) fol- 
lowed the initial separation for positive D/L. Thus, there is a distinct possibility that 
two totally different phenomena were being observed in the tests of reference 18, with the 
results at negative drag-lift ratios being a merely local phenomenon resulting from the 
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interaction of the model field with the poorly conditioned boundary layer that would be 
presumed to exist at such low tunnel speeds. 

The disagreement at negative D/L between the present theoretical study and the 
results of reference 18 probably cannot be resolved until additional experimental evidence 
becomes available. On the other hand, it is certain that any disturbances resulting in 
gross deformations of the wake obtained at CLhb = 3 and large negative D/L must 
occur far downstream of the model; consequently, any effects on the measured model 
performance should be minor. 

WIND-TUNNEL TESTS AT EXTREME CONDITIONS 
Performance in Free Air 

As is evident in figures 3 to 30 and 36 to 125, the flow near a model is decidedly 
distorted in the wind tunnel (compared to free air) when the wake deflection is sufficiently 
great. Even without flow reversal on the floor, small distortions (refs. 26 and 27) can 
have significant effects. The gross distortions of the wake shown by experiment (refs. 5, 
19, and 20) and indicated by the present theoretical treatment when the direction of wake 
flow is reversed, can, and generally will, produce such large effects on the measured 
wind-tunnel data that even gross trends in the data must be viewed with skepticism. 

Tests under these conditions are pointless and should be avoided. If data are required 
in this regime, which generally encompasses only the very low-speed end of transition 
flight, several alternatives appear to be available. These include tests of a smaller 
model in a larger tunnel, the use of a tunnel with more favorable proportions (fig. 34), or 
outdoor testing with one of the free-flight or test-bed techniques. 

The point at which wind-tunnel tests of a particular model should be abandoned can 
be determined at least roughly from the parameter xj|b given in figure 34. Although 
reference 18 appears to yield more restrictive limits at large negative D/L than fig- 
ure 34, it should be noted that reference 18 is not concerned with the effect on force mea- 
surements. The low wake deflections involved in the area of disagreement should result 
in placing the disturbances too far behind the model to have any significant effect on force 
measurements. 

Wake distortion and its effect on the overall flow in the tunnel is not the only crite- 
rion affecting V/STOL wind-tunnel data. Wall corrections, such as those of references 2 
to 4, can also become very large and very nonuniform over the model thus effectively 
distorting the aerodynamic shape of the model (ref. 24) until it no longer represents the 
geometric model for which data are desired. Indeed, since the flow both far forward of, 
and far behind, the model must be parallel to the walls, a large average correction always 
implies a large degree of nonuniformity (ref. 6). For large models, the errors associated 
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with nonuniformity of wall interference may be more restrictive than the onset of extreme 
wake distortion. 

The fact that the wall-interference theory (refs. 2 to 4) actually predicts the cor- 
rect effects on the overall flow confirms its general applicability. Numerous experi- 
mental studies corroborate this result. Indications to the contrary, presented in refer- 
ence 33, are disposed of in reference 6. 

The applicability of the theory of references 2 to 4 is not limited, as suggested in 
reference 18, to the initial appearance of flow reversal, or separation, on the floor, since 
the linearly inclined wake of the theory represents a wake which stagnates on the floor at 
least as well as it represents a wake which curves rearward without touching the floor; 
indeed the theory even predicts the flow reversal. The initial rollup of the front of the 
wake on the floor may not have much effect at the model since the image of the rolled-up 
wake in the floor will tend to cancel the effect of the real wake distortion at the model 
location. Instead, the limit of applicability of the theory is the point at which the conse- 
quent strength and position of the rolled-up vorticity have become such as to affect the 
data significantly. This point is, of course, that found experimentally by Rae (refs. 5 
and 20). 


Performance in Ground Effect 

The limiting conditions presented in figure 34 do not necessarily apply if the intent 
of the wind-tunnel tests is to obtain performance in ground effect. Examination of fig- 
ures 10 to 30 leads to the observation that the flow in the tunnel will differ little from that 
in ground effect provided that the model span is modest in comparison with the tunnel 
width. The wake distortions and rollup will still occur, of course, but these distortions 
are present in ground effect also. Some other limiting condition will eventually be 
reached when the rolled-up vorticity is essentially at the model station when it starts up 
the wall; however, the requisite conditions will be substantially more severe than those 
of figure 34. 

If the model is reasonably near the floor, interference at the model for the 
remaining three boundaries becomes small in comparison with the interference caused 
by the floor alone. (See ref. 2 and also table n for figs. 36 to 65.) Thus, neither the 
average wall effects nor their distribution will normally limit the extent of ground- effect 
tests. Neglect of the effect of the remaining walls is inherent in the data reduction for 
many tests and has been shown to yield good correlation with flight tests even for 
extraordinarily severe conditions (e.g., ref. 22). 

At short distances from the floor it is necessary to provide some means of elimi- 
nating the boundary layer on the floor since this boundary layer and its effect on the flow 
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field are not representative of true flight over a stationary surface. One successful 
means of treating the boundary layer is discussed in references 31 and 32. 

Flow Surveys 

Flow surveying, whether quantitatively or qualitatively (flow visualization), under 
severe conditions presents special problems since the desired information may be that at 
a point well removed from the model itself. Thus, figure 34, which applies only to force 
measurements at the model itself, is not applicable for determining limiting conditions 
for flow surveys. Even for mild conditions, such as in figure 3, the flow angles far down 
the wake are substantially different from those in free air. Under more severe condi- 
tions, such as in figures 7 to 9, there may be little resemblance between the flow in the 
tunnel and the flow in free air. Flow studies in this range of conditions (ref. 16) are use- 
ful only for examining wall effects per se. 

On the other hand, the flow in ground effect can be surveyed successfully to a far 
greater degree in the tunnel. As before, it is necessary for the model to be relatively 
small in span with respect to the tunnel and for the model relatively near the floor. In 
addition, care must be taken with the boundary layer on the floor. 

Under moderate conditions, flow surveys can be conducted quite successfully for 
V/STOL models. For example, calculations for the rotors of references 9 and 11 indi- 
cate that negligible wall interference was present. This result was obtained only because 
of the modest rotor size, the modest wake deflection, and the limited region surveyed. 

CONCLUDIN G REMARKS 

This study, based upon the theory of NASA TR R-124 and R-302 and an extended 
version of the theory of NASA TR R-71, indicates that the theory predicts the correct 
nature of the wake distortions which have been found to occur when testing V/STOL 
models under extreme conditions. The forward portion of the wake, which for moderate 
conditions flows rearward with the rest of the wake, will actually turn and flow forward 
if the overall wake deflection is sufficiently great. The forward-going portion of the wake 
will then roll up into a standing vortex on the floor, and its subsequent behavior appears 
to be that required to obtain the limiting conditions for V/STOL testing as found by Rae 
(Journal of Aircraft, May- June 1967). 

Theoretically, the predominant factor determining the initial tendency of the wake 
to flow forward appears to be the wake skew angle. Considerations of wake rollup indi- 
cate that model size and height should also be involved. A theoretical dependence on 
angle of attack (or, more generally, drag-lift ratio) is indicated; however, the theoretical 
result might be somewhat altered if wake curvature were considered. 
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Flow conditions such as to induce ceiling separation in a closed tunnel, and extraor- 
dinary distortion of the free streamlines that are the real open boundaries in a closed- 
on-bottom-only tunnel, are also indicated by the theory. In all cases, however, flow 
reversal on the floor is indicated as occurring prior to these latter effects. 

There is considerable doubt that wind-tunnel data obtained for conditions beyond 
Rae's limits are valid for free-air conditions. Valid ground-effect data can be obtained 
well beyond Rae’s limits provided certain conditions of model size and height above the 
floor are met. Some treatment of the boundary layer on the floor will be necessary as 
well. 

The effect of the calculated flow distortions is so great as to cast doubt on the 
effectiveness of wind-tunnel flow studies at the low-speed end of the transition range 
unless the results are intended to represent only the flow in ground effect. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 18, 1970. 
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APPENDIX A 


DERIVATION OF OPERATING PARAMETERS IN TERMS OF 
THE WAKE SKEW ANGLE 

Symbols 

A aspect ratio, b^/S 

Am momentum area of lifting system 

b span of lifting system 

Cjj induced-drag coefficient, D/qS 

Cl lift coefficient, L/qS 

CLhb lift coefficient based on tunnel cross-sectional area below and within model 

span, L/qhb 

D induced drag 

h height of model above floor 

L lift 

n ratio of induced velocities in final wake to induced velocities at lifting system 

1 9 

q dynamic pressure, — pV 

dt 

S wing, or rotor disk, area 

Uq mean, or momentum theory, value of longitudinal induced velocity (in 

X-direction) at lifting system, positive rearward 

V forward velocity 

Vr absolute value of the resultant velocity at lifting system 
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APPENDIX A — Continued 

Wp mean, or momentum theory, value of vertical induced velocity (in 

Z-direction) at lifting system, positive upward 

Wjj reference velocity, that value of Wq which would be obtained in hovering 

with D = 0 

X,Y,Z Cartesian coordinate system fixed in model, X positive rearward, Z posi- 
tive vertically, and Y to form a right-hand system 

a angle of attack, positive nose-up 

p mass density of fluid medium 

X wake skew angle, the angle from the negative Z-axis to the centerline of the 

wake at the model, positive rearward 


Theory 

Derivation of V/wp.- From figure 129, it may be seen that 


V + uq _ y u 0 
~Wp - " Wp ■ w 0 


(Al) 


On the assumption that the total mass flow affected by the aerodynamic system is 
equivalent to that which would pass through area A m with uniform velocity 


Vr = |V + up) 2 + (-w 0 ) 2 
forward thrust) are 


the lift and induced drag (which may be negative to represent 


L = pA m V R (-nw 0 J 

(A2) 

D = pA m V R (-nup) 

(A3) 


Dividing equation (A3) by equation (A2) yields 


so that equation (Al) becomes 


^0 _ D 
w 0 L 


tan x = 


V_ 

Wp 


D 

L 


(A4) 


(A5) 
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or 


V_ 

w 0 


= -(tan X+ 5 
L 


(A6) 


In the case of a rotor, the resultant force vector will be, in general, essentially 
perpendicular to the tip-path plane; thus, for a rotor 

— = tan a (A7) 

L 

Substitution of equation (A7) in equation (A6) yields 


V_ 

w 0 


-(tan x + tan o') 


(A8) 


Derivation of and Cj).- Define a reference velocity w^ as that vertical 

induced velocity which the system would require to hover with zero speed and induced 
drag and the same values of n and A m that pertain in forward flight. From equa- 
tion (A2) ^since w^ is equal to Wq under the aforementioned conditions^ 

L = pA m |w h |^-nw h j (A9) 


Solving for w^ yields 


w h = 



(A10) 


where the minus sign is required, since in the present coordinate system, a negative 
induced velocity (downward) is required to produce a positive lift (upward). 


For wings and rotors, n takes the value 2 (that is, in the far wake, the induced 
velocities are twice those at the lifting system) and the momentum area A m is essen- 
tially that of a circle circumscribing the lateral tips (A m = 7 rb^/ 4 ). Using these values 
and the definition of lift coefficient = L/qsV equation (A10) may be rewritten as 


w h 

V 


IC L (p/2)v 2 S 
2pl^. v 2 



(All) 
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Now the aspect ratio A of the lifting system is 



so that equation (All) becomes 



From equation (14) of reference 34 


| X | = co„ 



or 


Wq 

w h 


\Jcos x 


Thus, from equations (A13) and (A 15) 


X = .[ ffA 

w 0 w h w Q '^C L cos x 


Now equate equations (A6) and (A 16) to obtain 


(A12) 


(A13) 


(A 14) 


(A 15) 


(A 16) 


7tA I. D 

_ = - tan v + — 

C L cos x V L 


(A17) 


The solution of equation (A17) for is 


C L = 


7tA 


tan X + cos x 


(A 18) 


For a rotor, which has an aspect ratio of 4/7r, equation (A18) becomes 


C 


L = 


(tan x + tan a) cos x 


(A 19) 
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And, again for a rotor, in consequence of equation (A4) and the definitions of the 
coefficients 


C D = C L tan a = 4 — (A20) 

(tan x + tan a;) z cos x 

Derivation of Cl^.- Reference 18 defines limiting minimum test speeds in terms 

of a lift coefficient C Lhb based upon that portion of the cross-sectional area of the tun- 
nel contained beneath the wing span, thus, from equation (A18) 



APPENDIX B 


DERIVATION OF ROTOR INDUCED VELOCITIES IN 
WIND TUNNELS AND GROUND EFFECT 

Symbols 

The analysis presented in the following pages is valid only for an angle of attack of 
zero. In consequence, the wind-tunnel axes are parallel to axes in the rotor tip-path 
plane. Extrapolation to other angles of attack must be attended with great caution, since 
the direct extrapolation results in inconsistencies in x, z, a, wq, p, X, and x- 

B wind-tunnel semiwidth 

H wind-tunnel semiheight 

h height of rotor above floor 

m,N,n integers 

p integer, p = 1 for u and p = 0 for v and w 

q integer, q = 0 for closed tunnel and q = m + n for closed-on-bottom-only 

tunnel 

R rotor radius 

R c distance from point at (x,y,z) to edge of rotor disk at i//, 

||r 2 + + y2 + z2 - 2R(x cos xf/ + y sin 4/) 

r radius of vortex cylinder 

u,v,w components of induced velocity along the X-, Y-, Z-axes, positive when 

directed along positive direction of axis 

V forward velocity 

Wq induced velocity at center of uniformly loaded rotor (equal to mean, or 

momentum theory value of induced velocity normal to disk (ref. 35)), 
positive upward 
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X,Y,Z Cartesian axes centered in rotor, X positive rearward in plane of rotor, 

Z positive upward normal to plane of rotor, Y in plane of rotor to form 
a right-hand system 

x, y,z distances along X-, Y-, and Z-axes 

a angle of attack, positive nose-up (a = 0° throughout present analysis) 

y wind-tunnel width-height ratio, B/H 

£ semiheight of wind tunnel divided by height of model above floor, H/h 

X rotor inflow ratio, (V sin a + wq) j fiR 

p. rotor tip-speed ratio, V cos a/S2R 

v vectorial induced velocity, ui + vj, + wk 

a ratio of rotor radius to tunnel semiwidth, R/B 

X wake skew angle, angle between negative Z-axis and centerline of wake at 

rotor, positive rearward, tan _1 (-p/X) (for a=0°) 

\f/ rotor azimuth angle, positive counterclockwise from downstream position 

O rotor rotational speed 

Subscripts: 

G in ground effect 

T in tunnel 

t for triangular disk-load distribution 

trunc for truncated wake 
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Theory 

Uniform dis k- l oad distribution .- As noted in an earlier section, the wake of a uni- 
formly loaded rotor may be represented as a skewed cylindrical surface of uniformly 
distributed vorticity as in figure 130. Reference 10 obtains the induced velocities 
v(-m + vj + wk) of this wake as (note f(^) = 1 in ref. 10) 

' 2l R C \ 

1 (r + b- <= 03 xjcos X d X 

2tt 


u 

w 0 



R c R c / x\ 7 

-R R- + ( C 0 SX -Rj SinX + R C0SX 


R 

f + -^ cos xjsin X dx 


IT + fcos X - |)sin x + § cos x| 


R 


(Bl) 


(B2) 


P 

_w_ _ _ 1 _ l 

W 0 27 T \ 


2 77 " 


1-1 

fx y . s 

-C° s X + | sln X 

\ He . 

' + R~ Sin X cos X 

dx 

Rc 

r" 

r : + ( cosx -r) 

Sin X + R C0s X 



(B3) 


where equations (Bl) and (B2) have been corrected for an error in sign in equations (11) 
and (13) of reference 10.* 

Equations (Bl) to (B3) serve to define the rotor field in free air. It is possible to 
superpose the results to obtain the field in ground effect as indicated in figure 131. The 
first step is to translate the velocity field downward and rearward to its point of intersec- 
tion with the floor and then subtract the resulting field from the original free-air field to 
obtain the field of a truncated wake. Thus, if the free-air field is written symbolically as 


v _ v lx 2. 
w 0 WqIr’ R’ R .) 


the field of the truncated wake becomes 


(l\ =±k.z.^ 


tanx,£,t + fe 


\w Q j w n (R’ R’ R/ w n (R R — A ’ R’ R • R / 

7 trunc 


(B4) 


(B5) 


Incorrectly ascribed to sign convention in reference 6. 
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It is important to note that this truncated wake will coincide with the portion of the 
wake above the floor only if a = 0° . For any other value of a, the technique is inappro- 
priate and it would be necessary to return to reference 10 and perform the integrations 
along the wake to the floor, this distance being a function of azimuth angle and angle of 
attack. This effort does not seem warranted since the systems of references 2 to 4 will 
provide adequate results for most practical purposes. 

The next step in obtaining the flow field in ground effect is to rotate the field of the 
truncated vortex cylinder about the longitudinal axis of its intersection with the ground 
(fig. 131) and add the original and truncated fields. Thus, since the v and w compo- 
nents of the field are reversed in sign during the rotation, the total induced field in ground 
effect is 



■v fx y z\ v fx h , „ y z h \ 
w 0 (r- R’ Rj " w 0 (r " R U X ’ R’ R + Ry 




(B6) 


where p is an integer equal to 1 for the u component and is zero for the v and w 
components. For calculations with x = 90°, the wake never reaches the floor so that the 
second and fourth terms of equation (B6) may be neglected. 

Provided that the wind-tunnel floor is closed, the flow field in the wind tunnel is 
then obtained by reflecting the ground-effect field an infinite number of times in the hori- 
zontal and vertical boundaries, as in figure 132, to obtain 


/ \ oo oo T x 

4- £ ^ r ,q ^’ 

7 x m=-°° n=-°° v *- 


" H ' 2my | ‘ 4n 


cry 


cry 


-(-D* 4 



X 

- -i tan x 

I- - 2my * ■ 

- 4n + — 

V 

H 

e 

H H 


w 0 

_ 

cry ’ 

c ry J 

cry 



r~ 




V 

X 

JH 

- 2m y 
H r 

|-4n 

+ ! 

w 0 

oy’ 

cry 9 

cry 





V 

lx 

H " 

| tan x 

w 0 

oy ’ 


y 

H 


cry 


2 my - 4n + | 

? ~ 


cry 


(B7) 
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where q is an integer equal to m + n if the tunnel is closed on the bottom only, and q 
is equal to zero in the case of a closed tunnel. 

Notice that at x = 90°, the truncation of the wake occurs infinitely far behind the 
model. Numerical calculations at x = 90° can therefore simply delete the second and 
fourth terms of equation (B7). 

Equation (A6) (appendix A) yields a relationship between V and Wq which may be 
used to convert equations (B6) and (B7) to yield the field in terms of the free-stream 
velocity rather than the mean rotor- induced velocity; that is, for a = 0° 


_V_ 

w 0 


-tan x 


(B8) 


Thus equations (B6) and (B7) become 



(B9) 



(BIO) 


the 


If the total flow throughout the field is desired, it is only necessary to add ^ = 1 
u component of equations (B9) and (BIO). 


to 


Triangular disk-load distribution . - In the case of a triangular disk loading, the rotor 
is represented by a series of concentric skewed cylinders precisely as in reference 9. 

(See fig. 133.) The outermost cylinder, which has a radius equal to that of the rotor, must 
have a vorticity strength equal to 1.5 times the vorticity strength of a uniformly loaded 
rotor if the overall lifts of the rotors are to be equal (ref. 9). The 10 inner cylinders are 
evenly spaced in radius (starting with r/R = 0.05) and each is assigned a vorticity 
strength of -0.15 in order to reduce the vorticity in a stepwise manner to zero at the hub. 
Since the dimensions of the fields of the various cylinders differ only by constant factors 
depending upon the individual cylinder radii, the induced velocities for the entire system 
of cylinders representing a rotor with a triangular disk-load distribution may be written 
immediately as 




_L 

w 0 


/x_ I. 
r_ iL Ik 

r ’ _r ’ _r 
\R R R 


(Bll) 
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where 


r_ 

R 


N_ 

10 


0.05 


In order 

wind tunnel, it 
and (B7). 


to obtain the fields of the triangularly loaded rotor in ground effect or in a 

is only necessary to use (-^-) rather than in equations (B6) 

\ 0/t 0 


Conversion to tunnel centered co ordinates.- The equations in this appendix are all 
developed in terms of an axis system centered in the rotor, whereas in the main text it is 
more convenient to treat an axis system centered in the tunnel. The only difference is in 


the X-axis which is translated through a distance 



Thus, in terms of the tunnel 


coordinates the values for insertion in the equations of this appendix are 


x = x 


T 




y = y T 








(B12) 
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APPENDIX C 


FORTRAN PROGRAM FOR COMPUTING THE FLOW IN A WIND TUNNEL CONTAINING 

A LIFTING ROTOR 


THIS PROGRAM WAS WR ITT FN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE RFQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON T HE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 1? DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


INPUT WILL BE FOUND AT ADDRESS 1 (LINE (C19)I, IN FORMAT 100 (LINE 
(C 1071). THE REFERENCE ORIGIN IS THE CENTER OF THE ROTOR, WHICH IS RE- 
QUIRED TO BE AT AN ANGLE OF ATTACK OF ZERO. OUTPUT LOCATIONS (X, Y, ANO 
ZT) ARE WITH RESP C CT TO A REFERENCE ORIGIN SHIFTED VERTICALLY TO THE CENTER 
OF THE WINO TUNNFL . THE REQUIRED INPUT VARIABLES ARE: 


LI LOAD INDICATOR? LI=1 FOR UNIFORM DISK-LOAD DISTRIBUTION, L I =2 

FOR TRIANGULAR OISK LOAD DISTRIBUTION. (THE USE OF LI=2 
IS NOT RECOMMENDED BECAUSE OF EXCESSIVE RUNNING TIME. I 

ZETA SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT OF ROTOR ABOVE FLCOR . 

GAMMA WIDTH-HEIGHT PATIO OF WIND TUNNEL. 

SIGMA RATIO OF ROTOR DIAMETER TO WIND-TUNNEL WIDTH. 

CHI WAKF SKFW ANGLE, DEG. 

LOCATE DISTANCF OF CALCULATION PLANE TO RIGHT OF ORIGIN, BEHIND ORIGIN, 
OR ABOVE FLOOR (DEPENDING ON MAPTYPE), NONDI MENS ION AL I ZED 
WITH RESPECT TO THE TUNNEL SEMIHEIGHT. 

MAPTYPE INSERT LONGITUDINAL, LATERAL, OR FLOOR ( LEF T- JUST I F I F D » DEPEND- 
ING UPON DESIRED ORIENTATION OF CALCULATION PLANE. 


THIS PROGRAM IS NOT SUITABLE FOR CASES WHERE SIGMA IS ZERO SINCE NUMEROUS 
DIVISIONS BY ZERO WOULD BE REQUIRED IN SUBROUTINE <ROTUN>. CASES WITH CHI EQUAL 
TO EITHER ZERO OR 90-DEGREES ARE REJECTED. THESE CASES WOULD REQUIRE A DIVISION 
BY EITHER ZERO OR INFINITY AT LINE (C B5). 
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PROGRAM ROTWT ( INPUT , OUTPUT ,TAPE5=I NPUT, TAPE6=0UTPUT , PUNCH ) 

C* * *** *** *** *** ***** *************** ******************** ***************** 
C* * 

C* THIS MAIN PROGRAM SETS UP DO-LOOPS TO SCAN THE FLOW IN THE * 

C* DESIRED RFGION OF THE WIND TUNNEL* IT OBTAINS THE INDUCED FIELD * 

C* IN TERMS OF the mpaN INDUCED VELOCITY OF THE ROTOR AS WELL AS THE * 

C* COMPLETE FLOW IN TERMS OF THE FREE-STREAM VELOCITY. A PUNCHED * 

C* CARD OUTPUT OF THSE LATTER VELOCITIES IS PROVIDED FOR SUBSE- * 

C* QUENT PROCESSING TO PRODUCE VECTOR PLOTS OF THE FLOW. (SEE * 

C* APPENDIX D.) RFSULTS ARE OBTAINED SIMULTANEOUSLY FOR FREE AIR, * 

C* GROUND EFFECT, THE CLOSED WIND TUNNEL, AND A TUNNEL CLOSED ONLY * 

C* ON thp BOTTOM. ALL THRFE COMPONENTS OF VELOCITY ARE OBTAINED. * 

C* * 

C** **************************** ************** ** ** ***# ** ** ** *** *** ****** * 
COMMON LI , ZETA, GAMMA, SIGMA, X,Y,Z ,UOW , VOW , WOW , VELO( h , 3 > 

DIMENSION VEL(6) ,V(3) 

REAL LOCATf 

DATA < VEL ( II , 1 = 1,61 /4HU / WO* 4HV/WC ,4HW/W0 ,4HU/ V ,4HV/V ,&HW/V / 

1 READ (5,100) LI ,ZETA, GAMMA, SIGMA, CHI, L OC AT E , M APTYPE 

IF (EOF,?) 900,2 

2 IF (LI. EG. 2) GO TO 3 
ALO AD-10H UNIFORM 

GO TO 4 

3 AL.OAD=l°HTR T ANGUL AR 

4 WRITE (6,101) S TGMA, GAMMA, ZETA, ALOAO, CHI 

PUNCH 100, LT , ZFTA, GAMMA, SIGMA, CHI , LOC ATE , M AP T Y PE 
IE f MAPTYPE.EQ.*HL ATE) WRITE (6,103) LOCATE 

IP (MAPTYPF.EQ.4HL0NG) WRITE (6,104) LOCATE 

IE (MAPTYPE.EO.^HFLOO) WRITE (6,108) LOCATE 

WRITP (6,10?) 

I QU I T - 2 1 

IF (MAPTYPE.F0.4HFinO) IQUIT=20 
DO 11 1 1=1, 1 QUIT 

I c (MAPTYPE.EO.^HFLOO) GO TO 12 
ZTs-1 .1+C.1*FLQAT( I 1) 

Z= ZT + 1 . -1 . / Z FT A 

IE (MAPTYPE.EO.^HLATE) GO TO 12 

Y = L OC AT E 
TEND = 25 
GO TO 13 

12 X = LOC A T E 
TEND=L?.*GAMMA+1.0 

IF ( MAPTYPE .EO.&HL ATE) GO TO 13 

ZT=-1.0+LOCATE 

Z=-l .O/ZETA+LOCATP 

XF=TAN(CH 1*0.01 7 4532 92? 199) /ZETA 

I X F= X F 

ISTART=-?+TXF 

X=( FLOAT! I STAR T ) *0. 2)+C . 2*EL0AT(I1) 

13 DO 11 T ? =7 , TEND 

IF ( MAPTYPF.NE. ^HLONG) GO TO 15 
X=-2.2+0 .2*FL0AT (!? J 
GO TO 16 

15 Y=-O.l+D.1*FLOAT( 12) 

IF (12. EG. TEND) Y = G AMM A 

16 CALL ROTUN (CHI ) 

WRITE (6,106) X , Y , Z T 


(C 1) 
(C 2) 
(C 3) 
(C 4) 
(C 5) 
(C 6 > 
(C t> 
(C 8) 
(C 9) 
(C 10) 
(C 11) 
(C 12) 
(C 13) 
(C 1*H 
(C 15) 
(C 16) 
(C 17) 
(C 18) 
(C 19) 
(C 20) 
(C 21) 
(C 22) 
(C 23) 
(C 2M 
(C 25) 
(C 26) 
(C 27) 
(C 28) 
(C 29) 
(C 30) 
(C 31) 
(C 32) 
(C 33) 
(C 34) 
(C 35) 
(C 36) 
(C 37) 
(C 38) 
(C 39) 
(C 40) 
(C 41) 
(C 42) 
(C 43) 
(C 44) 
(C 45) 
(C 46) 
(C 47) 
(C 481 
(C 49) 
(C 50) 
(C 51) 
(C 52) 
(C 53) 
(C 5M 
(C 55) 
(C 56) 
(C 57) 
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C ** *** ********* ** *** * ** *** * ** * * ****** ** ********** ****** ** ** **❖ *** *** ** ** 

c* * 

C* THE FOLLOWING INDUCED VELOCITIES ARE IN TERMS OF WO: * 

C* * 

C* ********************************************************************** 
WRITE (6,102) <VELU),(VEL0(I,J),I=l,4>tJ = l,3) 

C** ********************************************************************* 


c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


NOTE THAT: 


1=1 

1=3 


IS 

IS 


IN 

IN 


FREE AIR 
CLOSED TUNNEL 


1=2 

1=4 


IS 

IS 


IN 

IN 


GROUND 

TUNNEL 


EFFECT 

CLOSED 


ONLY ON BOTTOM 


J=1 , 2 , 3 ARE UfV, AND W COMPONENTS OF INDUCED VELOCITY 


* 
* 
* 
* 
* 
* 
* 
* 

C,** ********************************************************************* 
IF (CHI.EO.PO.'M GO TO 11 
I c ICHI.EQ. 0 . 0 ) GO TO 11 
DO 5 1=1 ,3 

5 V(I) = VELU+ 3 ) 

C** ********************************************************************* 

C* * 

C* CONVERT TO PATIO WITH RESPECT TO V, NOTE V/WO=TAN( CHI ) * 

C* * 

C** ****************** ******** ******************************************* 
00 6 1 = 1 ,* 

DO 6 J = 1 , 3 

6 VFL 0 ( I , J)=-VFLPU , J)/TAN(CHI *0 . 0174532 Q 25199 ) 

C '.** ****** ** **** ****************** ***************** ********************** 

c* * 

C* ADD PR EE* STRE AM VELOCITY TO U COMPONENT * 

C* * 

C** ************ ********************************************************* 
00 7 1 = 1 , & 

7 VELOtl , 1 ) = 1 . 0 +VEL 0 ( I , 1 ) 

WRITE ( 6 , 102 ) ( V( J ) ,( VELO( I, J) , 1 = 1 , 4 ) , J=l, 3 ) 

d** ***************************** ***** *********************************** 

C* * 

C* PUNCH OUTPUT FOP SUBSEQUENT PROCESSING * 

C* * 

C* * **************************************************** ***************** 
PUNCH 10 -% ( X , Y, Z% ( (VELOCI ,J) , 1=1 , 4 ) ,J = 1 , 3 ) ) 

11 CONTINUE 
GO TO 1 

f ** **************************** ***************************************** 

c* * 

c* FORMATS * 

C* * 

C * * ** * *** * * * * * ** ** * * * * * * * ** * ** * *** * * * ** * * * * * * * * * *** ** * * * * ** * * * * * * * ** * * 

100 FORMAT ( II ,Fo # 3 t ^P! 0 . 3 , A 4 ) 

101 FORMAT ( 1 H 1 / / / /SO X* VEL OC I TI E S NEAR A LIFTING ROTOR*/// 

133 X* SIGMA =*fp. 3 t iOX*GAMMA = *F 8 . 3 , 10 X *Z ET A =*F 8 . 3 // 

237 X,A 1 %* DISK-LOAD D I S TR I BUTI 0 N *1 0 X*CH I =*F 8 . 3 //) 

102 FORMAT ( 5 X , A* • 4 F 2 5 . A ) 

103 FORMAT ( A 8 X*TR AN SVER SE SLICE AT X/H =*F 9 . 3 //) 

( 47 X*L 0 NGITUDINAL SLICE AT Y/H =*F 9 . 3 //) 

< 26 X*IN FREE AIR* 11 X*IN GROUND EF FECT* 8 X*I N CLOSED TUNNEL * 
C-O-B TUNNEL */ 26 X, ll(lH-) , 11 X, 16 ( 1 H- ) , 8 X , 1 6 ( 1 H- ) , 1 0 X , 15 ( 1 H- 


104 FORMAT 

105 FORMAT 
1 1 OX * I N 
2 ) ) 

106 FORMAT 
FORMAT 
FORMAT 
STOP 
END 


107 

108 
999 


(//2X*X/H =*F7.3,f X*Y/H = *F7. 3 , 5 X* Z /H =*F7.3/) 
I7F10.4/8F10.4) 

(41X*H0RI ZONTAL PLANE AT Z /H =*F9.3* ABOVE FLOOR*//) 


(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

<C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

(C 

<C 

(C 

<C 

(C 

(C 

(C 

(C 

(C 

(C 

( c 

(C 

< c 

(C 

(C 

(C 

(C 

(C 

(C 

tc 

( C 
IC 

tc 

( C 
(C 
(C 
(C 
(C 
(C 
(C 
(C 


58 ) 

59 ) 

60) 
61) 
62) 

63) 

64) 

65) 

66 ) 

67) 

68 ) 

69) 

70) 

71) 

72) 

73) 

74) 
75 ) 

76) 

77) 

78) 

79) 

80) 
81) 
82) 

83) 

84) 

85) 

86 ) 
87) 
88 ) 

89) 

90) 
91 ) 
92) 
93 ) 
94) 
95 ) 

96) 

97) 

98) 
90 ) 

100 ) 

101 ) 

102 ) 

103) 

104) 

105) 

106) 

107) 

108) 

109) 

110 ) 
111 ) 
112 ) 

113) 

114) 

115) 

116) 

117) 

118) 

119) 

120 ) 
121 ) 
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SUBROUTINE ROTUN (ANGLE) 


(C 

122) 

C* ******************************* ********************** ***************** ( £ 

123) 

c* 




* (C 

124) 

c* 


THIS SUBROUTINE OBTAINS THE INDUCED VELOCITIES INCLUDING 

THE * (C 

125) 

c* 


REAL ROTOR FOR THE IMAGE SYSTEMS REPRESENTING 

THE VARIOUS BOUNDARY* (C 

126) 

c* 


CONDITIONS. 


* (C 

127) 

c* 




* (C 

123) 

C** ******** *** ****** * * ********* ******** ****** *************************** ( C 

129) 



COMMON LI ,ZETA, GAMMA, SIGMA, X,Y,Z ,UOW, VOW, WOW, 

VELO ( 4, 3 ) 

(C 

130) 



DO A I =1 ,4 


(C 

131) 



DO 4 J =1 , 3 


(C 

132) 


4 

VELO( I, J)=0.0 


(C 

133) 



IF ( ANGLE • c O. °0 • 0) GO TO 1 


(C 

134) 



T AN C = T AN ( ANGLES .0174532925199) 


(C 

135) 


1 

AX = X $ AY = Y SA7 = Z 


(C 

136) 



DO 2 M=l,7 


(C 

137) 



DO 2 N = l,7 


(C 

138) 



AM=FL OAT ( M )-4. 


(C 

139) 



AN=FL04T(N)-A # 


( c 

140) 



X=AX/(SIGMA*GAMMA) 


(C 

141) 



Y=(AY/(SIGMA*GAMMA) ) -2. * AM/S I GM A 


(C 

142) 



Z= ( AZ-4 . * AN ) / ( SIGMA*GAMMA) 


(C 

143) 



CALL RTRLOAD (ANGLE) 


(C 

144) 



VELO ( 3,1)=VEL0(3,1)+U0W 


(C 

145) 



VELO ( 3 , ? )=VEL0(3,2)+V0W 


(C 

146) 



VELO ( 3,3) = VELO ( 3,3)+ WOW 


(C 

147) 



VEL0(4,1)=VEL0(4,1 ) + ( (-1 . ) ** ( M+N ) ) *UOW 


(C 

148) 



VELO (4, 2) =VEL0(4 ,2) + ( ( -1 . ) ** ( M+ N ) ) *VOW 


(C 

149) 



VEL0(4,3)=VEL0(4,3)+( (-1 . ) **( M+N ) )*WOW 


(C 

150) 



IF (M.NE.4.0R.N.NE.4) GO TO 3 


( c 

151) 



VELO ( 1, 1 )=VHLO( 2,1 )=UOW 


(C 

152) 



VELO( 1,2) =VE LO ( 2 ,2 ) = VOW 


( c 

153) 



VEL0(1,3)=VEL0( 2,3)= WOW 


( c 

154) 


3 

IF ( ANGLF.EQ.90.0) GO TO 5 


(C 

155) 



X=X-TANC/(SIGMA*GAMMA*ZETA) 


(C 

156 ) 



Z=Z+1 ./(SIGMA*GAMMA*ZET A) 


(C 

157) 



CALL RTRLOAD (ANGLE) 


(C 

158) 



VEL 0(3,1 ) =VE LO ( 3 ,1)-U0W 


(C 

159) 



VELO( 3,2) =VELO( 3, 2) -VOW 


(C 

160) 



VEL0(3,3)=VEL0(3,3)— WOW 


(C 

161) 



VEL 0 ( 4 , ! ) =VELO(&,l)-( (-1.) **(M+N) ) *UOW 


(C 

162) 



VELO(4,2)=VELO( 4, 2)-( ( - 1 . ) ** ( M+N ) ) *V0 W 


(C 

163) 



VEL0(4,3) =VEL0(4,3)-( (-1 . ) ** ( M+N ) ) *WOW 


(C 

164) 



IF (M.NE.4.0R.N.NE.M GO TO 5 


(C 

165) 



VELO ( 2 ,1 ) =VE LO ( 2,1) -UOW 


(C 

166) 



VELO(2»2)=VELO( 2, 2)* VOW 


(C 

167) 



VEL 0(2, 3 )=VEL0(2 ,3) -WOW 


( c 

168) 


5 

X= A X / ( SI GMA*GAMM A ) 


(C 

169) 



Y=-Y 


(C 

170) 



Z=( -AZ+A.*AN) /( SIGMA*GAMMA)-2. / (SIGMA*GAMMA*ZETA1 

(C 

171 ) 



CALL RTRLOAD (ANGLE) 


(C 

1^2) 



VELO ( 3,1)= VELO (3,1) +UOW 


( c 

173) 



VELO( 3,2 )=VEL0(?,2)-V0W 


(C 

174) 



VELO( 3,3) =VFLO( 3 , 3 ) ~ WOW 


(C 

175) 



VELO (4,l)=VELQ(A,l)+! (-!.)**( M+N )) *UOW 


(C 

176) 



VELO ( 4,2 )=VELO(4,2)-( ( -1 . ) ** ( M+N ) ) *VOW 


(C 

177) 



VELO(4,3)=VELO(^,3)-( ( - 1 • ) **( M+N) ) *W0 W 


(C 

178) 



IF (M.NE.4.0R.N.NE.4) GO TO 6 


(C 

179) 



VELO12,1)=VEL0( 2,1)+U0W 


(C 

180) 



VELO(2,2 )=VEL0(2 ,2) -VOW 


(C 

181) 
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APPENDIX C - Continued 


VEL0(2, 3) =VELO( 2, 3) -WOW 
6 IF ( ANGLE. F0.°0.0) GO TO 2 
X=X-TANC/( SIGMA*GAMMA*ZETA) 

Z=Z+1 ./ ( $ IGMA*GAMMA*ZET A) 

CALL RTRLOAD (ANGLE) 

VEL 0(3*1 )= VELO ( 3,1 )-UOW 
VFLO(3,2) =VEL0n,2)+V0W 
VFL0(3,3)=VEL0(3,3) +WOW 
VEL 0(4,1 )=VEL0(4,1 )-( (-1 . )**(M+NM*UOW 
VFLO(4,2)=VELO(* f 2)+( <-!.)**( M+N))*VOW 
VEL 0(4, 3 )=VEL0(4,3)+( (-1 . ) ** ( M+N ) ) *WOW 
IF (M.NE.4.0R.N.NE.4) GO TO 2 
VEL 0(2,1 ) =V FLO ( 2 » 1 1— UOW 
VELO( 2,2) =VEL0(2 ,2) + V0W 
VEL 0 ( 2 * 3 ) = VFLO ( 2,3 ) +WOW 
2 CONTINUE 

X= A X SY= AY SZ=AZ 

RETURN 

END 


(C 182) 
(C 183) 
(C 184) 
(C 185) 
(C 186) 
(C 137) 
(C 188) 
(C 189) 
fC 190) 
(C 191) 
(C 192) 
(C 193) 
(C 194) 
(C 195) 
(C 196) 
(C 197) 
(C 198) 
(C 199) 
(C 200) 


SUBROUTINE RTRLOAD ( ANGL ) 

C** # )$C # # 3$C £ * # ####* 3* 5$C 5}C # * * * * * 5ft * * *5?C # * jjt # & * % * J$C £ 5$C # S^C * J?C & * 3$C * afe * * # # # # * }}t * * # * * * 


C* * 
C* THIS SUBROUTINE PERFORMS 4 SUPERPOSITION ROUTINE TO OBTAIN * 
C* THE FIELD OF 4 TRIANGULARLY LOADEC ROTOR IF REQUIRED. THE * 
C* UNIFORMLY LOADED CASF FALLS DIRECTLY THROUGH TO SUBROUTINE * 
C* ROTVEL. * 
C* * 


C* * •+#*#♦***#*♦*##* &# **#**^!(t ** * #*** *##*>!'* ** =(!*** ****** * 

COMMON LT,ZETA, GAMMA , S I GM A , X , Y , Z , UOW , VOW , WOW , VELO ( 4 , 3 ) 

ACM I = ANGL SUOWT=VOWT=WOWT=C .0 $AX=X $ AY= Y $AZ = Z $NRING=11 
STR = 1 • 0 *IF (LI.EQ.l) NRI NG=1 $D0 2 I = 1,NRING 
IF (LI.EO.l ) GO TP 1 

AT = FLO AT ( T ) SR =0. 1 *A 1-0 .05 SX=AX/R $Y = AY/R $Z=AZ/R $STR=-9.15 
IF (T.NE.ll) GO TO 1 SX=AX *Y = AY SZ = AZ $STR=1.5 

1 CALL ROTVEL ( ANGL) 

UOWT=UOWT+UOW*STR SVOWT=VOWT+VOW*STR $ WOWT=WOWT+ WO W*STR 

2 CONTI NUF UUOW=UOWT $VOW=VOWT SWOW^WOWT 
RETURN SEND 


(C 201 ) 
(C 202) 
(C 203) 
(C 204) 
(C 205 ) 
(C 206) 
(C 207) 
(C 208 ) 
fC 209) 
(C 210) 
(C 211) 
(C 212) 
(C 213) 
(C 214) 
(C 215) 
(C 216) 

tc 2m 

(C 218) 
(C 219) 


SUBROUTINE ROTVEL (ACHI) (C 22 0) 

C* * (C 222) 

C* THIS SUBROUTINE OBTAINS THE INDUCED FIELD OF A UNIFORMLY * <C 223) 

C* LOADED ROTOR USING THE EQUATIONS OF NASA TN D-394. SEE NOTE * (C 224) 

C* FOLLOWING LINE fC 2^6 ) . * (C 225) 

C* * (C 226) 

C*********************************************************************** ( c 227) 
COMMON LI , ZET4, GAMMA, S IGMA, X,Y, Z,UOWt VOW , WOW , VELO ( 4, 3) (C 228) 

REAL INCRF «TTGHT=C.OOC01 $RT=X*X+Z*Z $ I F ( RT . GT . 2 . 0 ) T IGHT = T I GHT/R T (C 229) 
INCRE = 1.P * I NT = *• *ITEST = 2 $ITST=1 $R AD=0 .01745329251 99 $IT = 1 (C 230) 
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APPENDIX C - Concluded 


CHI=ACHI*RAD $TF (Y.LT.0.0.AN0.ACHI.EQ.9O.0I TT=2 * IF ( IT . EQ. 2 I Y=-Y 
IF ( ACHT .LE.90 .0 I GO TO 1 *ITE$T=1 *CHI =( 180.0-ACHI) *RAD $Z=-Z 

1 U=V=W=U1 =V1 =W1 =UOW=VOW=WOW=C.O $FAC=4.0 SPSI1=0.0 $IEND=360 
IF ( ACHI .EQ.90.0 ) GO TO 2 $RW=Y*Y+( X+Z*TAN( CHI )) **2-1 .0 

IF ( ABS(RWI .LT. 0.0001 .AND. Z.LS.O.O) GOTO 3 

2 IF ( ACHI. NE. 90.0. OR. Y.GT. 1.0) GO TO 7 

IF (ASSm. GT. 0.00001. OR. Y.EQ.O.C) GO TO 7 SRR=X*X+Y*Y 
IF (RR.GT.1.0. AND. X.LT. 0.01 GO TO 7 

IF (RR.GT ,1 .0. AND.X . GT .0.0 I GO TO 4 SPS 1 1= 1 80 .0- ASIN (YT *GO TO 7 

3 ZZ=X+Z*TAN(CHI1 *PSI 1 = ATAN2 ( Y , Z Z I tGQ TO 7 

4 DANGL=90 .0-ASIN(YI/RAD $1 ANGL=DANGL t INCRE =D ANGL /FL 0 A T( I ANGL! 

INT=1 SFAC=1.0 $IEND=2*T ANGL+1 $G0 TO 7 

5 INT=2 «GO TO 7 

6 I NT =3 $IENO=363-4-*I ANGL $ I NCPE= (360 .0-*- . 0 *0 ANGL I / FL 0 ATI I ENO-1 ) 

7 00 19 1*1* 1 END tPSI = FLOAT( I)*R AO+PSI 1 

IF (INT.EO.l) °SI=(90.0— INCRE*FLOAT ( 1-1 ) ) *R AD 
IF (INT.E0.2) P$I = (9O.0+INCRE*FLOAT(I-in*RAO 
IF fINT.E0.3I RST=(90.0+2.0*rANGL+INCRE*FLOAT(I-l))*RAD 
SC = SIN(CHI ) $CC=COS( CHI ) $SP=SIN(PSI5 $CP=COS(PSI) 

RC = S ORT ( 1 ,C+X*X+Y*Y+Z*Z-2.0*( X*CP+Y*SP ) 1 
IF ( RC.LT. 0.00001) GO TO 9 $ 8=RC+ < CP-X ) *SC+Z *CC 
IF (ABStBI.l_T.TI GHT 1 GO TO 10*A=1.0-( X*CP+Y*SP)+RC*SC*CP 
C=CP*( Z*RC*CC) SD=SP*<Z+RC*CC) 

U1=C / ( 8*R Cl $V1=D/(B*RC) $W1=A/(B*RC) *G0 TO l 7 
9 U1 =V1 =WI =0.0 1 1 F ( ACHI .NE.90 .0 I GO TO IF SGO TO 17 

10 IF ( ACHI. EO. 0.0) GO TO 14 *IF ( ABS ( Y) . L T.O .00001 ) GO TO 13 

11 UB= V B =WH=0 .0 *00 12 J=l,2 *PSII = PS 1-0 . 1 *R AD 
IF (J.E0.2) PSI I=PST+0.1*R AD 

RC= SORT ( 1 ,0+X*X+Y*Y + Z*Z-2.0*( X*COS f P ST I) +Y* SIN ( P SI I ) ) I 
A=1 .0-( X*COS( PST T)+Y*SIN(PSI I ) ) +PC*COS(PSI I ) *SC 
B=RC+(COS(PSII )-X)*SC+Z*CC *C=COSf PSI I )*( Z+RC*CC) 

0=SIN(PSTI)*( Z+RC*CC ) $UA=C/(8*PC) $VA=D/(8*RC) $WA = A/( B*RC ) 

UB = UB + UA *V8 =VB+ VA *W8 = WB+WA 

12 CONTINUE *U1 =U8/2 .0 iVl=VB/2.0 *Wl*WB/2.0 
IF (ITST.E0.2) GO TO 16 *G0 TO 17 

13 V1=0.0 *U1=X*SC*CC/(1.C-X*CP) *W1=1 .0+X*SC*CP/< RC*( X-RC*SC) ) 

GO TO 17 

14 Ul=-X/Z t Vl=—Y / Z *W1=1.0 *G0 TO I 7 

15 AX=X * AZ=Z *Z=Z-100. SX=X-Z*SC/CC *ITST=2 *G0 TO 11 

16 U1=U1/2.C «Vl=Vl/2.0 *Wl=Wl/2.0 $X=AX *Z=AZ *ITST=1 

17 IF (INT.NE.&. AND.I.EO.IEND) FAC=1.0 SU0W=U0W+U1*FAC*INCRE 
V0W=V0W+V1*F AC* INCRE $W0W=W0W+W1*FAC* INCRE 

IF (FAC. EQ. 4.0) GO TO 18 SFAC=4.C *G0 TO 19 

18 FAC* 2.0 

19 CONTINUE * I F (INT.EO.l) GO TO 5 $IF (TNT.E0.2) GO TO 6 

20 UOW* U0W/1080.0 SVOW* VOW/IOSO.O *W0W=W0W/1 080 . 0 

C************** ******************************************** ************* 


c* * 

c* NOTE THAT U/WO AND V/ WO HAVE BEEN CORRECTED FOR AN ERROR * 
C* IN SIGN IN EQUATIONS (11) AND (13) OF NASA TN D-39*. * 
C* * 


£** ********************************************************************* 
IF IITEST.EO.1) 1—1 SIF tIT.EQ.2) GO TO 21 $G0 TO 22 

21 Y=-Y *VOW=-VOW 

22 RETURM tENO 


( 

1 

( 


( 


(C 231) 
(C 232) 
(C 233) 
(C 234) 
(C 235) 
(C 236) 
(C 237) 
(C 238 1 
[C 239V 
(C 240 \ 
(C 241) 
1C 242) 
(C 243) 
C 244) 
:C 245) 
C 246) 
[C 247) 
C 248) 
C 249) 
C 250) 
C 251) 
C 252 ) 
C 253) 
C 254) 
C 255) 
C 256) 
C 257) 
C 258 ) 
C 259) 
C 260) 
C 261) 
C 262) 
C 263) 
C 264) 
C 265 ) 
C 266) 
C 267) 
C 268) 
C 269) 
C 270 V 
C 271) 
C 27Z'/ 
C 273) 
C 274) 
C 275) 

C 276 ) 
C 277) 

C 278 ) 

C 279) 

C 280) 
C 231) 

C 282) 

C 283) 

C 284) 

C 285) 
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APPENDIX D 


FORTRAN PROGRAM TO PREPARE CALCOMP TAPE 
FOR PLOTTING FLOW VECTORS 


THE FOLLOWING PROGRAM IS THE MOST GENERAL OF THOSE USED TO PREPARE THE FLOW 
VECTOR FIGURES PRESENTED IN THIS PAPER. SINCE THE SOFTWARE ASSOCIATED WITH THE 
PLOTTERS IN OTHFR COMPUTER INSTALLATIONS MAY VARY WIDELY, NO ATTEMPT IS MADE 
HEREIN TO PROVIDE THE LIBRARY SUBROUTINES FROM THE LANGLEY RESEARCH CENTER OF 
NASA. SINCE IT IS BELIEVED THAT MOST COMPUTING CENTERS WILL HAVE EQUIVALENT 
SUBROUTINES, THE FOLLOWING DESCRIPTIONS MERELY STATE THE FUNCTION OF THE SUB- 
ROUTINES AS USED HEREIN: 


AXES 

CALPLT 

CIRCLE 

LINE 

NOTATE 


DRAW A TIC-MARKED AXIS, WITH VALUES AT THE TIC-MARKS, AND 
IDENTIFY THE AXIS. 

MOVE PEN TO SPECIFIED X-Y COORDINATES. 

DRAW A SEMICIRCLE COUNTERCLOCKWISE FROM INITIAL X-Y COORD- 
INATES. 

DRAW A CONTINUOUS LINE THROUGH AN ARRAY OF POINTS. I TWO 
ADDITIONAL FINAL POINTS SPECIFYING SCALES ANO THE 
ORIGIN. MUST ALSO BE GIVEN.! 

DP AW ALPHANUMERIC INFORMATION FOR LABELLING. 


NUMBER CONVERT A FLOATING-POINT NUMBER TO ALPHANUMERIC INFORMA- 
TION, AND THEN DRAW FOR LABELLING. 


P ARROW DRAW AN ARROW FROM ONE SET OF X-Y COORDINATES TO ANOTHER 
WITH AN ARROWHEAD OF SPECIFIED LENGTH. 


INPUT TO THE PROGRAM WILL BE FOUND AT LINES «D 33), (D 35), AND (D 72). 
WITH THE EXCEPTION OF THE FIRST CARD, WHICH PROVIDES A FIGURE NUMBER, THE 
INPUT DATA IS PRECISELY THAT OBTAINED AS OUTPUT FROM THE PROGRAM OF APPENDIX C. 


PROGRAM CYL PLOT ( I NPUT , OUTPUT ,T APE5 = INPUT . T AP E6 = OUT PUT ) 

C* $$$$$$$ $$$$ $$$ ’Quf'X* $$$$ ** * $ 5(( $ $ 


C* * 
C* THIS PROGRAM PREPARES A CALCOMP TAPE FOR PLOTTING THE FLOW * 
C* VECTORS THROUGHOUT A PLANE NEAR A ROTOR IN AN ASSORTMENT OF WIND * 
C* TUNNEL CONFIGURATIONS. THE BASIC FLOW IS CALCULATED SEPARATELY * 
C* IN THE MANNER OF NASA TR R-71, BUT WITH THE LONGITUDINAL AND * 
C* LATERAL VELOCITIES ADDED BY THE USE OF THE EQUATIONS OF NASA TN * 
C* D-39A. PLANES MAY BE LONGITUDINAL, TRANSVERSE, OR THE FLOOR. * 
C* * 


C## ***#*#*# 

DIMENSION VELO(A-,3) ,XC(6) ,YC(6) , XD (75 ) , YD< 75 ) , S AV 1 ( 2 l , 25 ) , 

1 SAV2(21,25) ,SAV3< 21,25) , SAVA (21,25) ,SAV5< 21 ,25) ,S AV6 ( 21 , 25 ) , 

2 S AV 7 ( 21 ,25) ,SAV8( 21 ,25 )»SAV9(2l*25) ,SAV10(21»25) ,SAV11(21 ,25) , 


( D 

15 

ID 

25 

( D 

3) 

i D 

41 

(D 

55 

<D 

65 

(D 

75 

( 0 

8) 

< D 

9) 

( 0 

10) 

(0 

11) 

( D 

125 

ID 

13) 

(D 

14) 
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APPENDIX D — Continued 


3 SAV12(21,25J ,SAV13(21»25) »SAV14(21,25) , S AV15 ( 21 , 25 * 

REAL LOCATE 

£***************************************#******************************* 
C* * 

C* INITIALIZE PLOTTER. * 

C* * 

C*sfc *$ ********** ****************** ******** ************************ ******* 

CALL CALCOMP 
CALL LEROY 

CALL CALPLT (0.,0.,3) 

CALL C ALPLT(0 . , 5 . ,3 ) 

PI=3.141f 9265358979 

£******************************************** ***$*********************** 
C* * 

C* READ IN FIGURE NUMBER AND THE VARIABLES FOR WHICH THE * 

C* VECTORS ARE COMPUTED. * 

C* * 

C** ********************************************************************* 

8 READ (5,107) FIGURE 
IF (EOF, 5) 99°, 1 

1 READ (5,100) LI, ZETA, GAMMA, SIGMA, CHI .LOCATE, MAPTYPE 

IF (EOF, 5) 999,2 

2 IF (LI.E0.2) GO TO 3 
AL0AD=1DH UNIFORM 

GO TO 4 

3 AL0AD = 1DHTR I ANGULAR 

C* *************************************************************** ******* 

C* * 

C* WRITE OUT the CASE BEING PLOTTED * 

C* * 

C* ********************************************************* ************* 

4 WRITE (6,105) SIGMA .GAMMA, ZETA , ALOAD, CHI 

IF (MAPTYPE.EQ.4HLATE) WRITE (6,103) LOCATE 

IF (MAPTYPF.EO.tHLONG) WRITE (6,104) LOCATE 

IF (MAPTYPE.EO.^-HFLOO) WRITE (6,106) 

WRITE (6,108) FIGURE 

C* ********************************************************* ************* 

c* * 

c* INITIALIZE DO-LOOPS AND POSITION AND VELOCITY DATA. * 

C* * 

C** ************ ********************************************************* 

IOUI T =21 

IF ( MAP TYPE. EO. *-HFLOO) IQUIT=20 

I END=2 K 

IF ( MAP TYPE . NE . 4HL0NG ) IEND=10.*GAMMA+1 . 

DO 40 1=1,21 
00 AO J =1 ♦ 25 

AO SAV1 ( I, J)=SAV2l I, J)=SAV3< I , J)=SAV4( I , J)=S&V5( I ,J)=SAV6(I , J)= 

1SAV7CI » J ) =S A V8( T , J) = S AV9 ( I , J ) = S AV10 ( 1 1 J ) =S 4V 11 (I , J ) = S AV12 t 1 , J ) = 
2SAV13CI, J)=SAV14( I , J ) =S AV15 ( I , J ) =0. 0 
£******************** *******#*********************:(«*****#**** ****** ***** 

c* * 

C* READ IN AND ORGANIZE POSITION AND VELOCITY DATA. * 

C* * 

C** ********************************************************************* 
DO 7 I1=1,I0UTT 
DO 7 12=1 ♦ I END 

READ (5,101) CXfYtZTf MVELOdy J)«I = lt41f J = lt3) > 

IF ( EOF f 5 ) 999, 5 

5 S AVI ( II , 1 2 ) =X 
SAV2( II ♦ 1 2 ) =Y 


(D 15) 
(D 16) 
(D 17) 
(D 18) 
(D 19) 
(D 20) 
(D 21) 
(0 22 ) 
(D 23) 
CD 24) 
CD 25) 
(D 26) 
(D 2 7 ) 
(D 28) 
(D 29) 
(D 30 
(D 31) 
CD 32) 
(D 33) 
(D 34) 
(D 35) 
CD 36) 
(D 37) 
(D 38) 
(D 39) 
(D 40) 
(D 41) 
(D 42) 
(D 43) 
(D 44) 
(0 45) 
(D 46) 
(D 47) 
(D 48) 
(D 49) 
(D 50) 
(D 51) 
(D 52) 
(D 53) 
(D 54) 
(D 55) 
(D 56) 

( D 5 7 ) 
(D 58) 
(D 59) 
(D 60) 
(D 61) 
CD 62) 
CD 63) 
(D 64) 
(D 65) 
(D 66) 
(D 67) 
(D 68) 
(D 69) 
CD 70 
(D 71) 
(0 72) 
(D 73) 
CD 74) 
CD 75) 
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APPENDIX D - Continued 


SAV3(I1*I2)=ZT 

6 SAV4(Il,I21= VELOI 1,1) 

S AV5 (11*121= VELOI 1 » 2 I 
SAV6(IL,I2)= VEL 0(1 ,3 ) 

SAV7(IL,T21= VEL0(2,1> 

SAV8U1 .121= VELOC2,2» 

SA V9 (11.121= VEL 0(2,3) 

SAV10CI1 * 12 ) =VEL 0(3 , 1 ) 

SAV1KI1.T2) =VELO( 3*2) 

SAV12tIl,I2)=VEL0(3,3) 

SAV1 3( n,I2»=VEL0(4,l ) 

SAV14(I1. , T2»=VEL0(*, 2» 

SAVL5( U , 121 =VEL0<^,?» 

7 CONTI NiJP 

IF ( MAPTYPE. EQ. 4HL0NG) GO TO 30 
IF (NAPTYPF.EO.AHFLOO) GO TO 150 

C** ********************************************************************* 
C* * 

C* THIS SECTION PLOTS A TRANSVERSE PLANE. * 

C* * 

C** **<i*At***AAA«i(r*«!(!A*A*!!.*AAAi|i***t»**ti|i»AAA**A**A***>|iA*ili*Ai(r*A*A***A**** 
SC ALE = 0.5 

DO 1000 I PL 0T = ! ,4 

GO TO (31,32,33,34), IPLOT 

31 CALL NOTATE (1.0,-2.25,0.14,154(4).- FREE AIR. ,0.0, 15) 

GO TO 35 

32 CALL NOTATE ( 1 .0 , -2 . 25 , 0 . 1 4 , 20H ( B ) .- GROUND EFFEC T . , 0 .0 , 20 ) 

GO TO 35 

33 CALL NOTATE ( 1 . 0 , -2 . 25, 0 . 14 , 20H ( C ) .- CLOSED TUNNE L . , 0 . 0 , 20 ) 

60 TO 35 

34 CALL NOTATE (0 . 5 , - 2 . 25 , 0 . 14 , 35H( D ) .- CLOSED-ON-BOTTOM-ONLY TUNNEL. 

1 ,0.0,35) 

35 GO TO (250,251,251,251), IPLOT 

250 CALL NOTATE (-1 . 5 ,-2 . 75 , 0 . 1 0 ,7HF IGURE ,0.0,7) 

CALL NUMBER (99°. , 90°. , 0.10, FIG UR E, 0.0,0) 

GO TO 252 

251 CALL NOTATE (1 .0,-2.75,0.14 ,7HPI GURE ,0.0,7) 

CALL NUMBER (°°° . , 9°° . , 0 . 14 , F I GURE , 0 .0 , 0 ) 

252 GO TO ( 36,37,37,33) , IPLOT 

36 CALL NOTATE (99° . , ° 99 . , 0 . 10 , 46H- FLOW VECTORS IN A TRANSVERSE PLAN 
IE AT X/H = , 0 . 0 , 46 ) 

CALL NUMBER (99° . , °°° . , 0 . 1 0 , L OC ATE , 0 .0 , 3 ) 

CALL NOTATE ( 999. , °9° . , C. 1 0 , 3 6H , CALCULATED USING VORTEX CYLINCERS 

1 ., 0.0, 36) 

CALL NOT A TP (-1. B,-3. 0,0. 10 , 6 3H LOCATION OF ROTOR AND INTERSECTION 
1 OF WAKE AND PLANE ARE S HOWN . , 0 .0 , 63 ) 

GO TO 30 

37 CALL NOTATE ( 999 . , 999 . , 0 . 14 , 1 2H- CONT TNUED . , 0 . 0, 1 2 ) 

GO TO 30 

38 CALL NOTATF ( 9°9 . , 999 . , 0 . 14 , 12H- CONCLUDED . ,0. 0 ,1 2) 

30 IF (IPLOT. NE.l) GO TO 253 

CALL NOTATE ( - 1 . 5 , -7 . 25 , 0 . 1 0 , 7HZETA = ,0.0,7) 

C** ****** ***************** ** ***** ***** 44*44 ***************** ****** 

C* * 

C* THIS NOTATION GROUP IS USED FOP ALL PLOTS. * 

C* * 

C** ********************************************************************* 
129 CALL NUMBFR ( °99. , °99. , 0. 1 0 , ZET A ,0. 0 , 3 ) 

CALL NOTATP (99°. ,9°9., 0.10,24, ,0.0,2) 

CALL NOTATE ( 99°. ,99°. , 0. 10 , 30 , 0 .0 ,-l ) 

CALL NOTATE (99° . , 99° . , 0 . 10 , 3H = ,0.0,3) 
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APPENDIX D — Continued 


CALL NUMBER ( 999. , 9«° 0 .1 0 , GAMMA ,0 . 0 , 1 ) 

CALL NOTATE (999 . , 999 . , 0 . 10 , 2H , ,0.0,2) 

CALL NOTATE (999 . ,999 . , 0 . 1 0 ,40 , 0 .0 , -1 ) 

CALL NOTATE ( 999 . , 999 . , 0. 1 0 , 3H = ,0.0,3) 

CALL NUMBER (999 . , 999 . , 0 . 10 , S I GM A ,0 .0 , 2 ) 

CALL NOTATE ( 999. , 99«. , 0. 10 , 2H , ,0.0,2) 

CALL NOTATE (999 . , 999 . , 0 .10 , 28 , 0 .0 ,- 1 ) 

CALL NOTATE ( 999. , °99. , 0. 10 , 1 2H = 0.0 DEG, ,0.0,12) 

CALL NOTATE (99° . , 999 . , 0 . 10 ,43 , 0 .0 ,-l ) 

CALL NOTATE ( 999. , 999. , 0. 10 ,34 = ,0.0,3) 

CALL NUMBER ( 999 . , 99° . , C . 10 , CHI , C .0 , 1 ) 

CALL NOTATE ( 999 . , 999 . , 0. 1 0 ,4H 0EG,0.0,4) 

CALL NOTATE ( 9°9 . , °99 . , 0 . 10 , 1 8H , UNIFORM L 0 ADI NG . , 0. 0 , 1 8 ) 

IF ( MAPTYPE. EO. 4HL0NG) GO TO 130 
IF (MAPTYPE. EO.^HFLOO) GO TO 1E6 
253 XST4RT=(GAMMA-1 .0)/SCAlE 

IF (GAMMA. GT. 1. 1 ) X S T AR T= ( G AMM A-2. 0 1 /SC AL E 
XSTART1 = XST ART— 0 .5 
AXI S=4. 

IF (GAMMA. GT. 1.1) AXIS=8. 

CALL AXES(XSTART ,4.5 ,0.0, AXIS,-1. ,.5 ,2.0,2.5,3HY/H, .20, 3) 

CALL AXES ( 4 • *- 1 • »0 • * 2 • » 0 . , 2 • ?» • 4 » 2. 5 * 1 5H VECTOR SCALE, V, .14,-15) 
CALL AXES(XSTAPT1 ,0. ,9C. ,4. ,-l. , .R, 2.0 ,2.5 , 3HZ/H, .20, 3) 

c* * 

C* this SECTION PLOTS THE ROTOR AND THE INTERSECTION OF THE * 

C* WAKE AND THE PLOTTING PLANE. * 

C * * 

C## #$$$$$#***$*$#$*$$$$*** ## # * 4=4* * j{c j{e # Jfie* # ## * * 

DO 64 JKL MN = 1 f 

HH=( 1.0-ZET A ) / ? ET A 

EDGE =- SI GMA*GAMMA 

IF (LOC ATE. LE. EDGE) GO TO 62 

DO 61 J J =1 1 73 

AJ = FL 0 AT C J J 1 — 1 .0 

XD ( JJ)=SIGMA*GAMMA*COS( AJ*PI/36. ) 

YD ( J J ) = ( SIGMA*GAMMA*SIN( AJ*PI/36. )- LOC ATE) / T AN \ CH I * 0 . 01 745 32 92 51 99 
1 ) +HH 

IF (YD(JJ).GT.HH) YD(JJ)=HH 

IF (YD(JJ).LT. (-1,0) .AND. I PLOT .NE.l) YD(JJ)=-1.0 
IF ( YD( J J ).LT. (-2.0) ) YD ( J J ) =—2.0 

61 CONTINUE 

XD ( 74 ) =-GAMM& 

YD ( 74 ) =- 1 .0 

XD(T5) = YD ( 75 ) =SCALE 

CALL LINE(XD,YD*73»1, 0,0,0) 

62 XD(1I =SI GMA#G AMM A 
XDC 2)=-SIGMA*GAMMA 
YD ( 1 I =YD ( 2 ) = HH 

XD ( 3 ) =-G A MM A 
YDf 3 ) =-l .0 
XD ( 4 ) = YD ( 4 ) = SC AL E 
64 CALL LINE ( XO* YD , 2 « 1 , 0 , 0 , 0 ) 

c * 4 * 4 e 4 < 4 c 4 f 4 c 4 c 4 c 4 c 4 c 4 c 4 c 4 c 5 ?C 3 ^*j*: 35 c A 4c# 4c 3 ?c ^ 4 c afc 35 c ##*# # * # £ # 5 fc * sf: # #: # * $ * # # # # # :# # # # sjc # # * j}c #: 5 ^e # # # # s*t # # # * 

c* * 

C* THIS SECTION PLOTS THE FLOW VECTORS. * 

C* * 

C * 4c 4«4c4e4c4c4c 4«4c4<4c 4c4e ###### * #### # ## X#****#*^# A 4c A* 4c4c#4c4c * >}c ###:#:# #### #4c ### ### ##### 

63 DO 70 M=l,21 
DO 70 N=1 , IEND 

GO TO (52,53,54,55) , IPLOT 
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APPENDIX D - Continued 


52 V*SAV5(M,N) 

W=SAV6( M,N) 

GO TO 56 

53 V=SAV3(M,N) 

W=S AV9(M, N) 

GO TO 56 

54 V=S AV11 ( M , N ) 

W=S A VI 2 ( M , N) 

GO TO *6 

55 V=SAV14( M,N) 

W-S AVI 5 ( M ,N ) 

56 Q=SQRT( V* V+W*W) 

ALNGTH = . 0“V$C AL E 

IF(Q.LT.".5) ALNGTH=( .7*q*.2USCal e 
XA= (GAMMA+SAV 2 (M ,N) ) /SCALE 
XB=(GAMMA+SAV2 (M,N)+0 .2 *V US CAL E 
YA=( 1 .0+SAV3(M,N) ) /SCALE 
YB=(1 .D+SAV3 (M,N)+0.?*WU SCALE 
CALL PARROWf X A , Y A , XB , YB ,-l ,ALNGTF) 

IF (N.EQ.l) GO TO 70 

C 5 !'^ ************ #***### * ** **************************** *************** 


c* * 

C* NOTE that SYMMETRY IS USED TO OBTAIN THE LEFT-HAND HALF OF * 
C* THE PLANE. * 
C* * 


^** ************ ********************************** ********** **#***,}(£*: #,<(#* 


X A= ( GAMMA-S AV2 ( M , N ) )/ SC AL E 
XS=(GAMMA-SAV2f H,N!-0. ?*V) /SCALE 
CALL HARROW (X A, Y A,X 8 , YB,-1 , ALNGTH ) 

7 0 CONTIN'F 

CALL CALPLT (30. ,0., -3) 

ICO'’' CONTINUE 
GO TO B 

C************************** ****** ********************* ***** ** ****** ***** 

c* * 

C* this SFCTION PRODUCES LONGITUDINAL PLANES. * 

C* * 

(l* ************* ************** ***************************************,^* cJ/ . if 

80 DO 2000 TPL0T=1,4 

GO TO (21 1 22 ,23 ,2M , IPLOT 

21 CALL NOTATE (*. 0 ,-2. 25 ,0. 1* , 15H < A) FREE AIR. ,0.0,15) 

GO TO 25 

22 CALL NOTATF ( *. 0 2 . 2 5 , C . 1* , 20H ( B > . - GROUND EFF E CT . , 0 .0 , 20 ) 

GO TO 25 

23 CALL NOTAJc t 4. 0 ,-? . 25 , C . 14 , 20H I C ) . - CLOSED TUNNEL 0 .0 , 20 ) 

GO TO 25 

2 A- CALL NOTATE < 3 . c , -2 . 2 5 , 0 . 1 ^ , 35 H < C ) .- CL OS E D-ON- BOTTO M-ONLY TUNNEL. 
1 ,0.0,35) 

25 GO TO (260,261,261,261), IPLOT 

260 CALL NOTATF ( 2. 0 ,- 2 . 75 , 0. 10 , ^HFI GURE ,0.0,7) 

CALL NUMBFR (9QP. ,999. ,0.10, FIGURE, 0.0,0) 

GO TO 262 

261 CALL NOTAT<= ( 4 . " f _? # - 75 , 0 . 14 , T HF IGURE ,0.0,7) 

CALL NUMBER ( oeo . , OOQ . ,0.1^, FIGURE ,0.0*0) 

262 GO TO (26,2^,27,28), IPLOT 

26 CALL NOTATF ( <=?9 . ,coc . , c . 1 C ,49H- FLOW VECTORS IN A LONGITUDINAL PL 
1 AN E AT Y/H = ,0.0, * 8 ) 

CALL NUMBFR (°99. , 9°9 . ,0*1C,LOCATE,0.0,3) 

CALL NOTATF (9°9.,999.,0.10,36H, CALCULATED USING VORTEX CYLINDERS 

1. ,0.0,36) 

CALL NOTATE ( 2 . 0 ,-3. 00, 0. 10 ,^9HRCT0R AND FRONT AND REAR EDGES OF W 
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APPENDIX D - Continued 


27 

28 
2Q 


130 


1AKE ARE SHOWN. *0 .0*49 I ID 259) 

GO TO 29 ID 260) 

CALL NOTATE (999 . , 999 . , 0 . 14 , 1 2H- CONT INUEO . ,0 .0 , 12 ) CD 261) 

GO TO 29 (D 262) 

CALL NOTATE ( 99° . ,999 . , 0. 14 , 1 2H- CONCLUDED ., 0. 0, 12 ) ID 263) 

IF 1 1 PLOT . NE.l ) GO TO 130 ID 264) 

CALL NOTATE ( 2 . 0 ,-3 . 25 , 0 . 1 0 , 7HZ ETA = , 0 . 0 , 7 ) ID 265) 

GO TO 129 ID 266) 

CALL AXEStO. ,4 .5 ,0 . 0 , 10 . ,-2 . , . 5 ,2.0 ,2 .5 , 3H X/H , . 20 , 3) ID 267) 

CALL AXES!. 5,-1. ,0 . , 2 . , 0 . , 2 .5 , .4,2.5, 15HVECT0R SCALE, V, .14,-15) ID 266) 

CALL AXES(-.5,D.,90.,4.C,-1.,.5,2.0,2.5,3HZ/H,.20,3) 10 269) 

SC ALE=D . 5 ID 270) 

C*4 fr******************************************************************** ( o 271 ) 
c* * (D 272) 

C* THIS SECTION PLOTS THE ROTOR AND THE WAKE EDGES. * <D 273) 

C* * (0 274) 

( Q 275 ) 

XCm = (TAN!CHI*0.0174532925199)/ZETA)-SIGM4*GAMMA ID 27fc ) 

YCtl) =YCI4) =-1.0 ID 277) 

IF (XC(1).LE.3.0) GO TO 41 ID 278) 

YCtl) =11. /ZETA)- (3.0 +SIGMA*GAMMA)/T ANtCHI *0.01745 329 25199)- 1.0 ID 279) 

XCt 1 )=3.0 ID 280) 

41 XC 1 2 ) =-S I GMA*GAMMA ID 281) 

YCf 2)=YC(3)=(1./ZETA)-1. ID 282) 

XC(3)=-XC(2) ID 283) 

XC(4)=(TAN(CHI*0. 917*532925199) /ZETA)+SIGHA*GAMHA ID 284) 

XCI5)=-2.0 (D 285) 

YC I 5 ) =- 1 . 0 ID 286) 

XC(6 ) =YC(6)=SCALE ID 287) 

IF ( XCI4) .LE.3.0) GO TO 42 ID 288) 

YC(4)= 11. /ZETA)- (3.D-SIGMA*GAMHA)/TAN(CHI*0.C174532925199)-1.C ID 289) 

XC ( 4) =3. 0 ID 290) 

42 CALL LINE(XC,YC,4,1,0,0,0) ID 291) 

C** fr********************************#***#************’)'****************** ( o 292) 
C* * (0 293) 

C* THIS SECTION PLOTS THE FLOW VECTORS. * (D 294) 

C* * ID 295) 



DO 99 N= 1 , 25 

<D 

297) 


DO 99 M=1 *21 

< D 

298) 


GO TO <81,82,63,84) IPLOT 

(D 

299) 

81 

U=S 4 V4 ( M , N) 

< D 

300) 


W=SAV6(M,N) 

< D 

301) 


GO TO «5 

<D 

302) 

82 

U=S A V7 ( M, N) 

( D 

303) 


W=SAV<MM,N) 

( D 

304) 


GO TO 85 

(D 

305) 

83 

U=SAV19(M,N) 

(D 

306) 


W=S A VI 2 ( M , N) 

(0 

307) 


GO TO 8* 

( D 

308) 

84 

U=S A VI 3 ( M , N) 

(0 

309) 


W=SA VI 5 ( M ,N ) 

( D 

310) 

85 

0= SORT< U*U+W#W ) 

( D 

311) 


ALNGTH=. 07/SCALE 

(D 

312) 


IF ( 0. LT .0.5) AL NGTH= ( .7 *Q* . 2 ) / S C AL E 

( 0 

313) 


XAM2.0+SAV1 (M,N) ) / SCALE 

< D 

314) 


YA=(1 .9+SAV3(M,N) ) /SCALE 

<D 

315) 


XB=< 2.0+ S AVI <M,N) +0.2 *U) /SCALE 

( D 

316) 


YB=(1.0+SAV3<M,N)+0.2*W)/SCALE 

<D 

317) 


CALL PARROW ( X A , Y A , X B , Y 8 ,-l , ALNGTH ) 

< D 

318) 

99 

CONTINUE 

<0 

319) 
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APPENDIX D - Continued 


CALL C M_ PIT (3*;. ,0. ,-3) < 

2000 CONTINUE C 

GO TO 8 f 

C ^ ************ **************** ***************** **************** ******** ( 
C* * ( 

C* THIS SECTION PRODUCES VECTOR PLOTS OF THE FLOW ON THE FLOOR. * ( 

C* * ( 

C** ** ****** *** ***** ******** ****** **** ************* ***** ********* ******** ( 

150 DO 3000 I °L0T=1 ♦ A. < 

GO TO ( 151,152,153,154) IPLOT ( 

151 CALL NOTATE f3.0,-l. 7 5,0.14,15H( A) .- FREE AIR. .0.0*151 ( 

GO TO 15* l 

152 CALL NOTATE (3.0»-l. 7 5»0.14»20H(B).- GROUND EFFECT • . 0 .0 . 20 1 < 

GO TO 155 ( 

153 CALL NOTATE < 3 . n , -1 . 7 5 , 0 . 1* , 20H( C ) .- CLOSED TUNNEL ., 0 .0 ♦ 20 ) ( 

GO TO 155 C 

154 CALL NOTATE l 2 . 5 , -1 . 7 5 f 0 . 14 , 35H < D) CLOSED-ON- BOTTOM-ONLY TUNNEL. < 

1 ?0 . 0 ,35 ) ( 

155 GO TO 1270, 271, 271, 2 7 1) , IPLOT < 

270 CALL NOTATE ( 1 .0 ,-2 . 25 , 0 . 10 ,7HF I GURE .0.0,7) ( 

CALL NUMBER f 9*9. , 999. , 0. 10 , F I G UP E t 0. 0 ,0 ) < 

GO TO 2 7 2 < 

2 7 1 CALL NOT A TF f 3 . C , - 2 • 25.0. 1 4 , 7 H F I G UR E ,0.0,7) ( 

CALL NUMBFR (99 C .,9 C 9., 0.14, FIGURE, 0.0,0) ( 

272 GO TO ( 171, 172, 1 7 2, 173) , IPLOT ( 

1 7 1 CALL NOTATE < 9°9 . , 99* . , 0 . 10 , 2 3H- FLOW VFCTORS IN PLANE OF FLOOR,, ( 

1 0.0,33) ( 

CALL NOTATE t 9°° . , ooq . , q . 10 , 36H CALCULATED USING VORTEX CYLINDERS ( 
i. ,0.0,36) t 

CALL NOT ATF < 1 .0 ,-2 . 50 , 0 . 10 , 50H I NTERSECTI ON OF WAKE WITH PLANE OF ( 
1 FLOOR IS SHOWN. ,0.0, c 0) < 

GO TO 1^ ( 

172 CALL NOTATE (9 =>9 . k99° . , 0 . 14 , 1 2H- CONT INUED . , 0 . 0 , 1 2 ) < 

GO TO 1 7 4 < 

173 CALL NOTATE (99° • ,°99 .,0 .1^,12 H— CONCLUDED ., 0 . 0 , 1 2) < 

174 IF fIPLOT.NE.il GO TO 156 < 

CALL NOTATE ( 1 . 0 , -2 . 7 5 , 0 . 10 , 7HZ ET A = ,0.0,7) ( 

GO TO 129 < 

156 IFLOOR=TAN(CHI*0 .01745 3 29 25 19 9 ) / ZET A ( 

I ST ART =- 2 + 1 FLOOR < 

STARTUP AT(ISTART) ( 

CALL AXFSfO. ,-.5 ,C.O, 8 . , ST ART , .5 , 2 .0 , 2 . 5 , 3 HX/H , . 20 ♦ - 3 ) ( 

CALL AX C S( 7 .,-1.,0.,2.,0.,2.5,.4,2.5,15HVECT0R SCALE, V, .14,-15) ( 

CALL AXPS (-.5 ,o. ,90. ,2.0,0. *.5, 2.0, 2.5, 3HY/H, .20, 3) < 

158 SCALE=.* ( 

C** ** * ******* ****** ****** ****** ********************************** ******* f 
C* * t 

C* PLOT TH*= INTERSECTION OF THE WAKE ON THE FLOOR. * ( 

C* * f 

C ************** ********** * ******* ********************************* ****** [ 

ROCIR=RPCIR=SIGMA*GAMMA/SCALE < 

XC I R= f (T AN(CHI*0. r *17^5329251 c 9)/ZETA)-START) /SCALE + POCIR < 

YCIR=0. < 

CALL CIRCLEf XCIR,YCIR,0 .0, 180.P,R0CIR ,RFCIR,3) < 
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APPENDIX D — Concluded 


C** frjjc#**:)!*#*** lit ** ********* sit***##*#** fr********************** **** ******* 

C* * 

C* THIS SECTION PLOTS THE FLOW VECTORS. * 

C* * 

Q * a$c *** #** # # £ * * * *#* # *** * * * ** ** * ******** * *** * ****** ** **** * * ****** ** ***« * 

DO 199 M=1 , 20 

DO 199 N = 1 , T END 

GO TO (160,161,162,163) IPLOT 

160 U=SAV4(M,N) 

V=S AV5 ( M , N ) 

GO TO 16* 

161 (J=S A V7 ( M , N ) 

V=S A V8 ( M, N) 

GO TO 16* 

162 U= S A V 1 0 ( M , N ) 

V=S A VI 1 ( M , N) 

GO TO 16* 

163 U=SAV13(M,N) 

V=S A V 14 ( M , N ) 

164 X= S A VI ( M , N ) 

Y=SAV2( M,N) 

G=SORT(U*U+V*V) 

ALNGTH=.0 7 /SCALE 

IF (0.LT.0.5) ALNGTW=(.7*Q*.2)/SCALE 

XA= ( X- ST ART) /SCALE 

XB= ( X-START+.2*U)/SCALE 

YA=( Y5 /SCALE 

YB = ( Y+0. 2*V ) /SC ALE 

CALL PAR ROW (XA,YA 9 XB,YB,-1, ALNGTH ) 

199 CONTINUE 

CALL CALPLT (2*.t0.*-3l 
3000 CONTINUE 
GO TO 9 

C****************************************************** ***************** 

c* * 

C* TERMINATE PLOTTING. * 

C* * 

£************************** ****************************£********:******** 

99Q CALL CALPLT(^.,^ 0 ,90o) 

£************************** ****************:***£**#****#$*#***A*:£*>i<*;9::$:*** 

c* * 

C* FORMATS. * 

C* * 

C# * *************** *********£**#*****##*#*###*##*#*#****#****££#*#*£**#** 

100 FORMAT (T1?F9.3»^F10.3, A4) 

101 FORMAT ( 7 F10 .4/8F10 • * ) 

103 FORMAT ( * 8X*TR AN S VER SF SLICF AT X/H =*F9.?//) 

104 FORMAT (47X*L0NGITUDTNAL SLICE AT Y/H =*F9.3//) 

105 FORMAT (1H1////F0X* VELOCITIES NEAR A LIFTING ROTOR*/// 

1 33X*S I GM A =*F8.3,10X*GAMMA =*F8 .3, 10X *7 FT A =*F8.3// 

237X,A10,* DISK-LOAD D I S TR I BUTT 0N*1 0 X*CH I =*F*.3//) 


106 

107 


FORMAT 

FORMAT 

END 


(54X*VEL0CITIES 

(F1C.3) 
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VECTOR SCALE » V 


CB J — GROUND EFFECT. 


Figure 3.- Flow vectors in the X-Z plane, calculated using vortex cylinders. The rotor and the edges of the wake are shown. 
£ = 0.769; y = 1.0; a = 0.60; a = 0.0°; x = 70.0°; uniform loading. 
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Figure 3.- Concluded. 
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C B J — GROUND EFFECT. 

Figure 4.- Flow vectors in the X-Z plane, calculated using vortex cylinders. The rotor and the edges of the wake are shown. 
£ = 0.769; n = 1.00; y = 1.000; a = 0.600; a = 0.0°; x = 60.000°; uniform loading. 
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Figure 4.- Concluded. 
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Flow vectors in the X-Z plane, calculated using vortex cylinders. The rotor and the edges of the wake are shown. 
£ - 0.769; n = 1-00; y = 1.000; o = 0.600; a = 0.0°; x = 20.000°; uniform loading. 
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(B).- GROUND EFFECT. 

Figure 9.- Flow vectors in the X-Z plane, calculated using vortex cylinders. The rotor and the edges of the wake are shown. 
Z, - 0.769; n = 1.00; y = 1.000; o = 0.600; a = 0.0°; x = 10.000°; uniform loading. 
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Figure 10.- Flow vectors in a transverse vertical plane at x/H = 3.571, 
intersection of the wake on the vertical plane are shown. £ = 0.769; 
loading. 
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Figure 10.- Concluded. 
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Figure 11.- Flow vectors in a transverse vertical plane at x/h 
intersection of the wake on the vertical plane are shown. £ 
loading. 
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2.250, calculated using vortex cylinders. The location of the rotor and the 
0.769; n = 1.00; y = 1.000; o = 0.600; a = 0.00; X = 60.000°; uniform 
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Figure 12.- Flow vectors in a transverse vertical plane at x/H - 1.550, calculated using vortex cylinders. The location of the rotor 
and the intersection of the wake on the vertical plane are shown. £ - 0.769; y = 1.0; o = 0.60; a = 0.0°; x = 50.0°; uniform 
loading. 
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1.091, calculated using vortex cylinders. The location of the rotor and the 
0.769; n = 1.00; r = 1.000; o = 0.600 ; a = 0.0°; X = 40.0000; uniform 
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Figure 13.- Concluded. 
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Figure 14.- Flow vectors in a transverse vertical plane at x/h = 0.750, calculated using vortex cylinders. The location of the rotor and the 
intersection of the wake on the vertical plane are shown. £ = 0.769; n = 1.00; y = 1.000; o = 0.600; a = 0.0°; x = 30.000°; uniform 
loading. 
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Figure 15.- Flow vectors in a transverse vertical plane at x/H = 0.464, calculated using vortex cylinders. The location of the rotor and the 
intersection of the wake on the vertical plane are shown. £ = 0.769; n = 1.00; y = 1.000; o = 0.600; a = 0.0°; x = 20.000°; uniform 
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. Figure 17.- Concluded. 
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Figure 18.- Flow vectors in the plane of the floor, calculated using vortex cylinders. The intersection of the wake 
and the plane is shown. C = 0.769; n = 1.00; y = 1.000; a = 0.600; a = 0.0°; x = 60.000°; uniform loading. 
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Figure 19.- Flow vectors in the plane of the floor, calculated using vortex cylinders. The intersection of the wake 
and the plane is shown. Z = 0.769; n = 1.00 : y = 1.000; a = 0.600; a = 0.0°; x = 50.000°; uniform loading. 
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Figure 19.- Concluded. 
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Figure 20.- Flow vectors in the plane of the floor, calculated using vortex cylinders. The intersection of the wake 
and the plane is shown. Z, - 0.76 9; n = 1.00; y = 1.000; a = 0.600; a = 0.0°; x = 40.000°; uniform loading. 
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Figure 21.- Concluded. 
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Figure 22.- Flow vectors in the plane of the floor, calculated using vortex cylinders. The intersection of the wake 
and the plane is shown. £ = 0.769; y = 1.0; a = 0.60; a = 0.0°; x = 20.0°; uniform loading. 
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Figure 24.- Continued. 
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Figure 28.- Isometric views of the flow vectors near a rotor, calculated using vortex cylinders. The rotor and the intersections of the wake 
with the X-Z plane and the plane of the floor are shown. Because of symmetry, only the vectors in the right-hand half of the plane of 
the floor are given. Z, = 0.769; y = 1.0; o = 0.60; a = 0.0°; x = 30.0°; uniform loading. 
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Figure 28.- Continued. 
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Figure 29.- Isometric views of the flow vectors near a rotor, calculated using vortex cylinders. The rotor and the intersections of the wake 
with the X-Z plane and the plane of the floor are shown. Because of symmetry, only the vectors in the right-hand half of the plane of 
the floor are given. ^ = 0.769; n = 1.0; r = L0; o = 0.60; a = 0.0° ; x = 20.0°; uniform loading. 
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Figure 20 - Continued. 
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Figure 30.- Isometric views of the flow vectors near a rotor, calculated using vortex cylinders.' The rotor and the intersections of the wake 
with the X-Z plane and the plane of the floor are shown. Because of symmetry, only the vectors in the right-hand half of the plane of 
the floor are given. C “ 0.769; n = 1.0; y = 1.0; o = 0.60; a = 0.0°; x = 10.0°; uniform loading. 
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Figure 30.- Continued. 
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Figure 30.- Concluded. 
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Figure 33.- Comparison of calculated flow field with tuft-grid study in a closed-on-bottom-only tunnel. Calculations are for 
plane of tuft-grid (0.13H to side of rotor). Calculated flow field assumes triangular disk-load distribution, y = 2.0; 
o = 0.25; x ~ 35°; Z = 1.0. 
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Figure 34.- Minimum-speed limits for V/STOL testing as derived from the data of reference 5. 
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Figure 36.- Concluded. 
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Figure 37.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 0.769; n = 1.000; y = 1.000; a = 0.600; a = 0.000°; x = 70.000°; uniform loading. 
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Figure 38.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. Z, = 0.769; n = 1.000; y = 1.000; o = 0.600; a = -20.000°; x = 70.000°; uniform loading. 
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Figure 39.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ - 0.769; n = 1.000; y - 1.000; o = 0.600; a = 20.000°; x = 50.000°; uniform loading. 
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Figure 39.- Concluded. 
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Figure 40.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. S = 0.769; n = 1.000; y = 1.000; o - 0.600; a = 0.000°; x = 50.000°; uniform loading. 
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Figure 42.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. Z, = 0.769; n = 1.000; y = 1.000; a = 0.600; a = 20.000°; x = 30.000°; uniform loading. 
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Figure 42.- Concluded. 
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VECTOR SCALE . V 


IB) — GROUND EFFECT- 

Figure 43.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 0.769; n = 1.000; y = 1.000; o = 0.600; a = 0.000°; x = 30.000°; uniform loading. 
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Figure 43.- Concluded. 
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Figure 44.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edqes of the 
wake are shown. Z = 0.769; n = 1.000; r = 1.000; o = 0.600; a = -20.000°; x = 30.000°; uniform loading. 
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Figure 45.- Flow vectors in the X-Z plane, 
wake are shown. £ = 0.769; n = 1.000; 


calculated using doublet strings. The rotor and the edges of the 
y = 1.000; o = 0.600; a - 0.000°; x = 10.000°; uniform loading. 
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Figure 45.- Concluded. 
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Figure 46.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 1.000; n = 1.000; y = 1.000; o = 0.600; a = 20.000°; x = 70.000°; uniform loading. 
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Figure 46.- Concluded. 
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(6).- GROUND EFFECT. 

Figure 47.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. Z, = 1.000; n - 1.000; y = 1.000; o = 0.600; a = 0.000°; x = 70.000°; uniform loading. 
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Figure 48. - Flow vectors in the X~Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 1.000; n = 1.000; y = 1.000; o = 0.600; a = -20.000°; x = 70.000°; uniform loading. 
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Figure 48.- Concluded. 
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(8).- GROUND EFFECT. 


Figure 49.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 1.000; n = 1.000; y = 1.000; o = 0.600; a = 20.000°; x = 50.000°; uniform loading. 
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Figure 50.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 1.000; n = 1.000; y = 1.000; a = 0.600; a = 0.000°; x = 50.000°; uniform loading. 
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Figure 51.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. C - 1.000; n - 1.000; y = 1.000; o = 0.600; a = -20.000°; X = 50.000°; uniform loading. 
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(8).- GROUND EFFECT. 


Figure 52.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. ^ = 1.000; n = 1.000; y = 1.000; o = 0.600; a = 20.000°; x = 30.000°; uniform loading. 
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( B ) .- GROUND EFFECT. 

Figure '53.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. C = 1.000; n = 1.000; y = 1.000; o = 0.600; a = 0.000°; x = 30.000°; uniform loading. 
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Figure 53.- Concluded. 
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IBJ.- GROUND EFFECT. 

Figure 54.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 1.000; n = 1.000; y = 1.000; a = 0.600; a = -20.000°; x = 30.000°; uniform loading. 
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Figure 54.- Concluded. 
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(B).- GROUND EFFECT. 


Figure 55.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 1.000; n = 1.000; y = 1.000; a = 0.600; a = 0.000°; x = 10.000°; uniform loading. 
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Figure 55.- Concluded. 
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f B) •- GROUND EFFECT. 

Figure 56.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 1.428; n = 1.0; y = 1.0; o = 0.600; a = 20.0°; x = 70.0°; uniform loading. 
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VECTOR SCALE, V 

(B).- GROUND EFFECT. 

Figure 57.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = 1.428; n = 1.0; y = 1.0; a = 0,600; a = 0.0°; x = 70.0°; uniform loading. 
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Figure 57.- Concluded. 
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Figure 58.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. Z = 1.428; n = 1.0; r = 1.0; o = 0.600; a = -20.0°; x = 70.0°; uniform loading. 
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Figure 59.- Flow vectors in the X~Z plane, calculated using doublet strings. The rotor and the edges of the 
wake are shown. £ = ‘1.428; n = 1.0; y = 1.0; a = 0.600; a = 20.0°; x = 50.0°; uniform loading. 
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Figure 59.- Concluded. 
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(B).- GROUND EFFECT. 

Figure 60.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.428; n = 1.0; y = 1.0; a = 0.600; a = 0.0°; x = 50.0°; uniform loading. 
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Figure 60.- Concluded. 
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Figure 61.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.428; n = 1.0; y = 1.0; a = 0.600; a = -20.0°; x = 50.0°; uniform loading. 
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Figure 61.- Concluded. 
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Figure 62- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, - 1.428; n = 1.0; y = 1.0; a = 0.600; a = 20.0°; x = 30.0°; uniform loading. 
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Figure 62.- Concluded. 
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(B).- GROUND EFFECT. 

Figure 63.- Flow vectors*in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z = 1.428; n = 1.0; y = 1.0; a = 0.600; a = 0.0°; x = 30.0°; uniform loading. 
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Figure 63.- Concluded. 
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(B)-- GROUND EFFECT- 


Figure 64,- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.428; p = 1.0; y = 1.0; o = 0.600; a = -20.0°; x = 30.0°; uniform loading. 
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Figure 64.- Concluded. 
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VECTOR SCALE . V 


(8).- GROUND EFFECT. 

Figure 65.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.428; n = 1.0; y = 1.0; a = 0.600; a = 0.0°; x = 10.0°; uniform loading. 
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Figure 65.- Concluded. 
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Figure 66.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.000; n = 1.000; y = 1.000; o = 0.600; a = 20.000°; x = 70.000°- tri- 
angular loading. 
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Figure 66.- Concluded. 
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Figure 67- Concluded. 
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Figure 68.- Concluded. 
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(B).- GROUND EFFECT. 

Figure 69.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, - 1.000; n = 1.000; y = 1.000; a = 0.600; a = 20.000°; x = 50.000°; tri- 
angular loading. 
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(B).- GROUND EFFECT. 


Figure 70.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.000; n = 1.0; y = 1.0; o = 0.600; a = 0.0°; x = 50.0°; triangular 
loading. 
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Figure 70.- Concluded. 
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(BK- GROUND EFFECT* 

Figure 71.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, - 1.000; n = 1.000; y = 1.000; 0 = 0.600; a = -20.000°; x = 50.000°; tri- 
angular loading. 
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(B)- GROUND EFFECT. 


Figure 72- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.000; y = 1.000; o = 0.600; a = 20.000°; x = 30.000°; tri- 
angular loading. 
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(D)- CLOSED-ON-BOTTOH-ONLY TUNNEL. 

Figure 72.- Concluded. 
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(B).- GROUND EFFECT. 


Figure 73.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edqes of 
the wake are shown, q = 1.000; n = 1.000; y - 1.000; o = 0.600- a = 0 000°- x = 30 000°- tri- 
angular loading. 
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Figure 73.- Concluded. 
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(B) - GROUND EFFECT. 


Figure 74.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.000; y = 1.000; o = 0.600; a = -20.000°; x = 30.000°; tri- 
angular loading. , 
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Figure 74.- Concluded. 
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Figure 75.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z = 1.000; n = 1.000; y = 1.000; o = 0.600; a = 0.000°- x = 10.000°- tri- 
angular loading. 
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Figure 76.- Concluded. 
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( B ) -- GROUND EFFECT. 


Figure 77.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.000; y = 1.000; a = 0.300; a = 0.000°; x = 70.000°; uniform 
loading. 
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Figure 78.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, - 1.000; n = 1.000; y = 1.000; o = 0.300; a = -20.000°; x “ 70.000°; uniform 
loading. 
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Figure 79.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 1.0; 0 = 0.300; a = 20.0°; x = 50.0°; uniform loading. 
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Figure 80.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1-000; T = 1-000; a = 0.300; a = 0.000°; x = 50.000°; uniform 
loading. 
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Figure 80.- Concluded. 
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VECTOR SCALE. V 


CB ) — GROUND EFFECT. 

Figure 81.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. C = 1.000; n = 1.0; y = 1.0; o = 0.300; a = -20.0°; x = 50.0°; uniform loading. 
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Figure 82.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.000; y = 1.000; a = 0.300; a = 20.000°; x = 30.000°; uniform 
loading. 
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( B) GROUND EFFECT. 


Figure 83.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.000; n = 1.000; y = 1.000; a = 0.300; a = 0.000°; x = 30.000°; uniform 
loading. 
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Figure 83.- Concluded. 
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IB) GROUND EFFECT. 


Figure 84.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.000; y = 1.000; o = 0.300; a = -20.000°; x = 30.000°; uniform 
loading. 
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Figure 84.- Concluded. 
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(B)- GROUND EFFECT. 


Figure 85.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.000; y = 1.000; a = 0.300; a = 0.000°; x = 10.000°; uniform 
loading. 
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C B 1 -- GROUND EFFECT. 

Figure 86.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 1.0; o = 0.900; a = 20.0°; x = 70.0°; uniform loading. 
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Figure 86.- Concluded. 
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Figure 87.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z = 1.000; n = 1.0; y = 1.0; o = 0.900; a = 0.0°; x = 70.0°; uniform loading. 
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( B ) — GROUND EFFECT. 

Figure 89.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z = 1.000; n = 1.0; y = 1.0; o = 0.900; a = 20.0°; x = 50.0°; uniform loading. 
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Figure 91.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z = 1.000; n = 1.0; y = 1.0; o = 0.900; a = -20.0°; x = 50.0°; uniform loading. 
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Figure 92.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z = 1.000; n = 1.000; y = 1.000; o = 0.900; a = 20.000°- x = 30 000°- uniform 
loading. 
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Figure 93.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000* n = 1.000; y = 1.000; a = 0.900; a = 0.000°; x = 30.000°; uniform 
loading. 
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Figure 93.- Concluded. 
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Figure 95.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.000; y = 1.000; a = 0.900; a = 0.000°; x = 10.000°; uniform 
loading. 
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Figure 99.- flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 2.0; o = 0.600; a = 20.0°; x ~ 50.0°; uniform loading. 
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Figure 99.- Concluded. 
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Figure 100.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 2.0 ; a = 0.600; a = 0.0°; x = 50.0°; uniform loading. 



266 



\ 0 



i ! 144-4 

VECTOR SCALE. V 


(CJ - CLOSED TUNNEL. 



A I 44-4-4 

VECTOR SCALE. V 


1 l 1 1 

D> O ^ — £> — > 
[>{>£> C> O O -O 
C> O £> 0-000 

{>£>[> 0^£> — £ — £► 
t> O D> E > O D> 4 

O {> t> t>— 1> O -O 

t> E> O t> t> 0 O 

E> O £> 0- 0 D> fo 
O O t> 0- 0 t> 4> 
{>£>{> 0-0 T> -O 

O O O 0-0 0- 0 

t> E> 0 O O O O 

O O O O -0-0-0 
O O O O l> t>C> 
!> O 0> O t> 0 t> 

O i> O O O D> O 

O O r> O D> i> O 
£>!>!> O l> t> O 
£>!>!>{> - O 0 0> 
[>[>£> O O 0 E> 
O O O O O 0 O 


CD) CLOSED-ON-BOTTOM-ONLY TUNNEL. 

Figure 100.- Concluded. 


267 



X/H 


-i L 



L 1 I I I 

0 2 3 4 5 

VECTOR SCALE , V 


"tA 

AA -"A. -"A ' 

''A ^A 
"A "A' 

"A "A ""A ~ 

'"A ''A ^A "A ^~A " 

"A ^ ^A 
''A "'A -"A ^A ^A~ 

'A "A "A ""A ~ 

''A ""A ~A ^A ~ 

^ A ^A "A 'A ~A~ 

^A ^A ^A ~ 

'A -A> A£.~ 


(AJ- FREE AIR- 


-l l 



1 -+ 4-1 ! I 

VECTOR SCALE * V 


(BA- GROUND EFFECT- 


Figure 101 - Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.000; n = 1.0; y = 2.0; o = 0.600; a = -20.0°; x = 50.0°; uniform loading. 
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Figure 102.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. ? = 1.000; n = 1.0; y = 2.0; o = 0.600; a = 20.0°; x = 30.0°; uniform loading. 
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Figure 103.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 2.0; a = 0.600; a = O.OP; x = 30.0°; uniform loading. 
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Figure 104.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.000; n = 1.0; y = 2.0; a = 0.600; a = -20.0°; x = 30.0°; uniform loading. 
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Figure 105.- Flow vectors in 
the wake are shown. C - 


the X-Z plane, calculated 
1.000; FI = 1.0; y = 2.0; 


using doublet strings. The rotor and the edges of 
o = 0.600; a = 0.0°; x = 10.0°; uniform loading. 
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Figure 106.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 2.0; a = 0.300; a = 20.0°; x = 70.0°; uniform loading. 
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Figure 107.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. ^ = 1.000; n = 1-0; y = 2.0; o = 0.300; a = 0.0°; x = 70.0°; uniform loading. 
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Figure 109.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 2.0; o = 0.300; a = 20.0°; x = 50.0°; uniform loading. 
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Figure 110.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 2.0; o = 0.300; a = 0.0°; x = 50.0°; uniform loading. 
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Figure 111.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. ^ = 1.000; n = 1.0; y '= 2.0; o = 0.300; a = -20.0°; x = 50.0°; uniform loading. 
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Figure 111.- Concluded. 
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Figure 112.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 2.0; a = 0.300; a = 20.0°; x = 30.0°; uniform loading. 
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Figure 113,- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. K = 1.000; n = 1.0; y = 2.0; a = 0.300; a = 0.0°; x = 30.0° ; uniform loading. 
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Figure 114.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. I = 1.000; n = 1.0; y = 2.0; o = 0.300; a = -20.0°; x = 30.0° ; uniform loading. 


294 





X/H 

-2 -1 0 1 2 3 4 5 

j i I I I | I I I I | ! I I I | t I I ! | I ! I I | T — T i I f J \ i 1 | 



VECTOR SCALE, V 


CC) CLOSED TUNNEL. 


X/H 

-2 -1 0 1 2 3 4 , 5 

j I I I I j ! i i i | i i i t j i i i l | i ! i i | i i 1 » j i i J i j 



VECTOR SCALE. V 


(DJ — CLOSED-QN-BOTTOrt-ONLY TUNNEL- 
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Figure 116.- Concluded. 
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Figure 117.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 0.5; o = 0.600; a = 0.0°; x = 70.0°; uniform loading. 
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Figure 118.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.000; n = 1.0; y = 0.5; o = 0.600; a = -20.0°; x = 70.0°; uniform loading. 
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Figure 119.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 0.5; o = 0.600; a = 20.0°; x = 50.0°; uniform loading. 
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Figure 120.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 0.5; a = 0.600; a = 0.0°; x = 50.0°; uniform loading. 
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Figure 122.- Flow vectors in the X-Z plane, calculated using doublet strings, 
the wake are shown. ? = 1.000; n = 1.0; y = 0.5; o - 0.600; a = 20.0°; 


The rotor and the edges of 
X = 30.0°; uniform loading. 
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Figure 123.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. Z, = 1.000; n = 1.0; y = 0.5; o = 0.600; a = 0.0°; x = 30.0°; uniform loading. 
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Figure 124.- Flow vectors in the X-Z plane, calculated using doublet strings. The rotor and the edges of 
the wake are shown. £ = 1.000; n = 1.0; y = 0.5; a = 0.600; a = -20.0°; x = 30.0°; uniform loading. 
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Figure 125.- Flow vectors in the X-Z 
the wake are shown. C ~ 1.000; n 


plane, calculated using doublet strings. 
= 1.0; y = 0.5; a = 0.600; a = 0.0°; 


The rotor and the edges of 
X = 10.0°; uniform loading. 
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Figure 126.- Concluded. 
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Figure 127.- Lines of constant wake skew angle x in terms of c Lhb com P ared with the values for a closed tunnel 

extrapolated from reference 18 . 









(a) Closed tunnel. 

Figure 132.- Central portion of reflection pattern to obtain induced velocities in tunnel. Solid lines indicate closed boundaries; 
broken lines indicate open boundaries. Sense of vorticity is indicated by signs, + indicates sense same as basic tunnel model 
and - indicates sense opposite to basic tunnel model. Real boundaries are darkened. 
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(a) Wake as a nest of vortex cylinders. (b) Load distribution. 

Figure 133.- Superposition technique used to obtain the induced field of a rotor having axisymmetric nonuniform disk-load distribution. 
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